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JIAIN ROTOR FREE WAKE GEOMETRY EFFECTS ON BLADE AIR LOADS 
AND RESPONSE FOR HELICOPTERS IN STEADY MANEUVERS 
VOLUME I - THEORETICAL FORMULATION AND ANALYSIS OP RESULTS* 

By S. Gene Sadler 

Rochester Applied Science Associates, Inc. 


SUMMARY 


Rotor wake geometries were predicted by a process similar to 
start-up of a rotor in a free stream. An array of trailed and shed 
vortices was generated with vortex strengths corresponding to step- 
wise uniform radial and azimuthal blade circulations. The array was 
limited to an arbitrary number of azimuthal steps behind each blade, 
and the remainder of the wake model consisted of an arbitrary number 
of trailed vortices behind each array. Vortex elements were allowed 
to be transported by the resultant of the free stream and vortex- 
induced velocities . Blade angles of attack and relative motion be- 
tween the rotor and its previously deposited wake were influenced 
by maneuvers. Wake geometries were computed, as were wake- induced 
velocity influence coefficients for use in the calculation of non- 
uniform downwash in the blade loads calculations . 

Blade loads and response computations assumed steady-state 
(periodic) conditions, and included maneuver effects, the effects 
of nonuniform inflow, nonlinear airfoil characteristics, gravity 
loading, gyroscopic loading, and the response of flexible blades to 
the applied loads. The resulting nonlinear equations were solved 
by an iterative process, and computed data was compared with mea- 
sured data, where available. The measured turn maneuvers were 
approximately steady maneuvers, the pull-ups were not steady, and 
no roll data was available. An alternate vortex-blade interaction 
model was used for one pull-up case. 

Vortex-blade interactions were clearly indicated in the calcu- 
lated wake-induced velocities. Air loads calculated with the use 
of static aerodynamic stall characteristics did not reflect many of' 
these vortex-blade interactions, and were particularly unresponsive 
to high frequency angle of attack changes in the static stall re- 
gion of angle of attack for an airfoil representation which included 
detailed Mach number effects in that region. Effects of both angle 
of attack rate and compressibility are required to predict air loads 
for such conditions. Computed blade response quantities were con- 
sistent with computed air loads and the normal modes used to repre- 
sent the blades. 


*VOLUME II - PROGRAM LISTINGS is contained in NASA CR-2111. 



INTRODUCTION 


Maneuver capabilities currently represent an undesirable limit 
on the operational envelope for most helicopters. These limitations 
are often a result of severe oscillatory air loads and helicopter 
response. There is an increasing need for an understanding of the 
causative physical phenomena, and there has recently been an in- 
creased effort to develop tools and to perform experimental and 
theoretical investigations to gain an understanding of these phe- 
nomena and to develop methods to overcome the existing limitations. 

Accurate prediction of the aerodynamic flow field in which the 
helicopter rotor operates is critically important in the calcula- 
tion of helicopter air loads, blade response, and vibration levels. 
There are indications that vortex-blade interactions are often 
dominant sources of high frequency air lo.ads of sufficient magni- 
tude that severe oscillatory response results. Such interactions 
have been discussed in references 1 through 4, among others. Meth- 
ods for predicting the rotor's aerodynamic flow field have recently 
undergone extensive refinement, primarily due to the need for im- 
proved predictions and the availability of high-speed, large- 
capacity digital computers. These refinements have resulted in 
more realistically detailed models for calculation of nonuniform 
wake-induced inflow. Some of these models are discussed in refer- 
ences 5 through 9 . Corresponding blade air load and blade response 
models have been developed, so that compatible blade loading 
and response models are used. The improved models for the wake- 
flow prediction have provided a good qualitative approximation to 
average measured wake flows (refs. 5, 10) and are expected to per- 
mit improved prediction of higher harmonic air loads and blade 
response . 

Reference 5 discussed some of the difficulties associated with 
the prediction of vortex- induced inflow, and described a model 
which was developed to predict wake geometry and associated wake- 
induced flow and wake-induced velocity influence coefficients. 

This model showed good agreement with the measured average wake 
flows reported in reference 10, but was restricted to steady level 
forward flight operating conditions. The model discussed in this 
report is an extension of that discussed in reference 5, and in- 
cludes effects of steady maneuvers on both wake geometry and blade 
load and response calculations . The computer program which has 
been implemented for this investigation is a modification of that 
of 'reference 5, and is operational for single rotor helicopters 
only. Wake geometries, wake-induced velocity influence coeffi- 
cients, blade air loads and blade response were computed. Blade 
response (motions, shears, and moments) were computed in terms of 
a finite set (up to 12) of normal modes. The loads/response cal- 
culations were implemented for fully articulated rotors, teetering 
rotors or for rotors with (assumed) rigid hubs. 
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The blade loads/response model assumes periodiGity of loads 
and response, so Is inherently limited to analysis of steady- 
state level and maneuver flight conditions, and cannot be extended 
for analysis of non-steady flight conditions. The flight condi- 
tions analyzed herein include comparisons with some experimental 
cases which were not steady maneuvers, but which were intended to 
provide a measure of the acceptability of the programs as pres- 
ently implemented. 

The wake model program uses a forward-time integration method 
of solution. Vorticity is defined in terms of straight vortex 
elements with uniform properties along their lengths . The present 
wake program treats only steady-state level and maneuver flight 
conditions but could be extended to predict wake geometries , wake- 
induced inflow, and approximate blade loads and response for non- 
steady flight conditions as well, such as gust loading, transient 
maneuvers, etc. The extension of the wake program to include cal- 
culations of loads and response would necessitate a set of sub- 
routines which would compute blade loads and response with a for- 
ward time integration method. It is anticipated that such exten- 
sions would increase the core storage requirements significantly, 
but would not greatly increase running times . 


ANALYSIS 


The computer programs were implemented as two main programs 
and preliminary calculations corresponding to convenient divisions 
in the theoretical representation of the helicopter rotor system and 
its wake, referred to hereafter as the model formulation and solu- 
tion of governing equations. The program usage flow is indicated 
by the diagram in figure 1. The blade loads and response program 
uses input generated by the wake program as well as normal modes 
and natural frequencies defined by a previously developed blade 
frequency program. Both the wake and the blade loads and response 
programs require input specifying operating conditions , trim con- 
trol settings, etc. These inputs may be obtained from experimental 
data or from appropriate preliminary calculations . The combined 
analyses and programs were directed toward prediction of blade 
air loads and response for helicopters in steady -maneuver flight 
conditions, in which the effects of a free wake and flexible blades 
have been included. (One flight condition calculation was per- 
formed with an empirically developed vortex-blade interaction model • 
implemented for use in defining aerodynamic lift and moment coeffi- 
cients . This was an exception to the usual procedure . Results of 
this case are noted at appropriate parts of the report.) As pres- 
ently implemented, only one rotor may be modeled, blade-wake and 
wake-wake interactions were allowed, but no mass nor elastic cou- 
pling through the hub between blades was allowed. (The computer 
programs used for forv/ard flight conditions could model two 
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arbitrarily located rotorSf but core storage requirements and other 
considerations resulted in the present Single-rotor model imple- 
mentation for maneuver conditions . ) All blades on the rotor were 
assumed to be identical, and all blades' motions, loads, etc., were 
also assumed to be the same. 


Wake Model and Formulation 

The wake geometry is calculated by carrying out a process 
similar to the start-up of a rotor in a free stream. The blades 
are located at a specified azimuthal and flapping position, without 
any wake. The blades then rotate through an azimuthal increment, 

A\jj , and shed and trail vortex elements of unknown strength, but with 
known positions. The strengths of the vortices that are shed imme- 
diately behind the blade are then determined, and include the effect 
of their own self-induced velocities. (If it is desired, an. esti- 
mate of the blade loads that result can then be determined, without 
the effects of blade flexibility being included.) All vortex ele- 
ment end points not attached to the blade are then allowed to trans- 
late as the blade is stepped forward for a time At with velocities 
as determined by the free stream and induced velocities. Here, 

At = Aip/Q, where Q is the rotational speed. This con5>letes the 
first step in the wake geometry calculation. 

Subsequent steps are similar. The blades are restated through 
an angle Ai|; to new (specified) flapping angles. The vortices at 
and just behind the blades have unknown strength, while those in 
the wake have known strength. The strengths of the vortices 
attached to the blade are then computed, induced velocities at 
vortex element end points are computed, and end points not attached 
to the blade are allowed to translate with the resultant velocity 
due to the free stream and induced velocities. (Again, if desired, 
approximate blade loads may be computed, once blade circulations 
are known.) After the blade circulation values are determined for 
a given azimuthal location of the blades, vortex- induced velocities 
are computed at all vortex element end points in the wake. The 
blades are then advanced through an azimuthal increment, Af, and 
each vortex end point which is not attached to the blades is allowed 
to translate for the time period A^/Q, with the resultant velocity 
due to the free stream and induced velocities and with the position 
changes due to maneuvers, if any. The entire computational process 
for the new rotor blade positions is then repeated. That is, new 
wake-induced and free stream velocities at the blades are computed 
and used to determine new blade circulations . Then new vortex end 
point locations are computed, etc. Each blade advance results in 
an additional set of shed and trailed vortices being added to the 
wake. The number of revolutions of wake retained for actual compu- 
tational purposes is specified by an input to the computer program. 
For practical use, the number of revolutions required is given by 

N^l/jJ IT . 
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In the manner described above, arrays of discrete shed and 
trailing vortices are generated iramediately behind the blades with 
strengths which correspond to approximate blade loads. These arrays 
have stepwise radial and asimuthal strength variations so that total 
circulation is conserved. The arrays of shed and trailing vortices 
which are generated immediately behind the blades are referred to as 
the full mesh wake. Comparisons of wake flows predicted using this 
wake model for the entire wake with experimental measurements in- 
dicated that retention of shed elements with a coarse mesh resulted 
in poor induced velocity predictions , and that use of a fine mesh 
increased running time to an unacceptable level. Therefore, the 
full mesh wake was used to represent the wake immediately behind 
the blades, and a modified wake model was developed and implemented 
for use in the representation of the remainder of the wake, as 
shown in figure 2. 

The modified wake consists of trailing vortices only, so vor- 
ticity is not conserved. The locations, circulations, core sizes, 
etc. of the trailing vortices in the modified wake are determined 
in terms of the full mesh wake's last sets of trailing vortex loca- 
tions , circulations , etc . The wake-induced velocities , wake distor- 
tions, and other calculations are essentially the same for both the 
full mesh wake and for the modified wake portions of the wake model. 

The wake and blade vortices are assumed to have finite-sized 
cores of rotational fluid, and the core sizes at the blade are con- 
trolled by an input parameter so as to be adjustable to improve 
agreement between calculated and experimental or theoretical core 
size data as desired. All elements are assumed to be straight for 
purposes of caleulations , except where induced velocity, calcula- 
tions are performed for an end point of a vortex element due to an 
element with one end at that point. At such a point the neighboring 
elements are assumed to be arcs with curvature determined by three 
appropriate end points . 

The basic mathematical formulations of the wake geometry pro- 
gram are those related to determination of blade motion, calcula- 
tion of vortex- induced flow, and calculation of blade circulations 
including self-induced effects. 

Blade motion. - The determination of approximate (steady and 
first harmonic) blade flapping motion is accomplished by using 
experimental data if available, or data from performance type cal- 
cuiations . If cyclic and collective pitch values and appropriate 
blade properties are known but blade coning and flapping motions 
are not known, a program option is available to compute approximate 
steady and first harmonic flapping motions internally. This option 
is described in detail in reference 5, and does not include maneuver 
effects . 
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Vortex-Induced flow . -Vortex induced flow is computed at blade 
aerodynamic load points for calculation of bound circulation, and 
at vortex element end points for calculation of wake distortions. 

The Biot-Savart law Csee ref. 11, for example) is the basic rela- 
■tionship used for ealculating vortex induced flow exterior to 
vortex elements' rigidly rotating cores. At calculation points 
located inside a rigidly rotating vortex, the vortex induced 
'Velocity is determined from the rigidly rotating core velocity at 
that point. Flow at a vortex element end point due to the adjoin- 
ing vortex elements accounts for curvature in the manner developed 
by Crimi in reference 7. At points where calculated velocities 
from an individual element are unrealistiGally high the velocities 
are restricted to be less than rotor tip speed; the total calculated 
velocity at a point is restricted to be less than a specified (by 
input) percentage of rotor tip speed. This type of limit is phys- 
ically realistic, since the application of a very high velocity to 
a calculation point for a finite time increment. At, results in a 
larger motion than V70uld result from a truly continuous type of 
calculation. That is, due to the finite nature of the model, cal- 
culation of wake distortions from vortex- induced velocities is 
unrealistic when very large velocities are permitted to remain over 
finite times, but are reasonable when velocities are limited. (In 
practice the limits were applied infrequently, as indicated by the 
program output.) The method of calculation of vortex- induced flow, 
i.e., the Biot-Savart law, is not directly influenced by the maneuver 
flight condition model requirements. Other effects of maneuvers on 
wake geometry, etc., are included in other parts of the model. 

Maneuver effects on wake position . -The wake-induced and free 
stream velocities, together with a velocity due to the relative 
rotation between the rotor and the wake which accompanies the 
steady-state, maneuvers are then computed to give each wake element 
end point its particular total velocity. The blades are then step- 
ped forward through an angle Arp, and the wake element end points 
are transported by their respective velocities for a time At, where 
At = Arp/^. Relative rotations and translations between the wake- 
oriented (inertial) coordinate system and the rotor-oriented 
coordinate system, as shown in figure 3, may be described by a 
coordinate transformation of the form 



where the first term represents the translation due to helicopter 
motion, the second term represents the rotation due to helicopter 
motion, and the third term represents the wake-induced transla- 
tions. Here (x, y, z) are the rectangular coordinate locations of 
the point. The definitions of R and S for the steady-state forward 
flight and maneuver flight conditions are given below for coordinated 
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turns ^ pull'-ups/ and rolls. 


^ For coordinated turns with, a rotation rate of S about a vertical 
axrs C+ for no se-^ to -right turns) 

(Rj . (A](PCA))(I>]R, 

|s| = (Aj j^P(A } j j-p Co} C2) 

where 

0 +COS 01 +sina 

s s 

10 0 
0 -sina cos a 
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-cosAsin^^ 
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c 
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sin$^ 
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c c 


cosA -1 
|Dj = sinA 
0 

|p(o)j = |P(A) j with A = 0 

R = turn radius , meters 
c 

= shaft tilt and positive aft, radians 
s 

= tan“^ (V|/R^g) 

= helicopter velocity, meters/second (see fig. 4) 

A = 6 At 

At = 2 tt/N r 
s 

N = number of azimuthal steps per revolution 

O 
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and otiier aymbols have been defined previously. For a steady roll 
about a roll axia at roll rate positive for right’- side -down 

roll, the transformation matrices are defined as 
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the roll center is located at in the shaft coordinate 

system, and 

= angle (about the y-axis) between the roll axis 
and the x~axis , radians 

= I At. 


Symmetrical pull-up at a rate with a radius R^ from the hub 
axis result in coordinate transformation matrices as defined by 
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sina^ 

0 R 

P 

1-cosa^ 


cosa^ -sina^ 

0 10 
sina^ 0 cosa.^ 

where = a^At. 

Steady climb or steady forward flight with climb rate w and/or 
forward velocity result in coordinate transformation matrices 
defined by 

^f 

0 

w 

c 


1 0 0 

S = 0 1 0 

[o 0 1 

The above forms of transformations are not limited by small angle 
approximations, so effects on wake location of maneuvers of any 
severity or time duration may be simulated. Effects of gravity 
loading and gyroscopic terms will be discussed in a later section, 
inasmuch as these maneuver effects primarily modify blade loads 
calculations . 
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Determination of bound circulation . -Blade circulations are 
dependent upon maneuver parameters, as determined by their aero- 
dynamic effects. The velocities V and U, normal and tangential to 
the plane normal to the shaft axis, respectively, are given approxi- 
mately by 


ttg) + w + r (,a^ cos ijj + 


cos sin i})) 


nr + sin ijj 


where r = radial position, meters, 

U = total tangential velocity, meters/second, 

V = normal velocity, meters/second, 
w = normal wake-induced velocity, meters/second, 

= pitch rate, positive nose up, radians/second, 
a„ - angle between roll axis and x-axis, radians, 

ri 

a = shaft axis angle with respect to the free stream, 

“ positive aft, radians, 

a = longitudinal first harmonic blade flapping, positive for 
flap up at ii)=0 , radians , 

= roll rate, positive right-hand side down, radians/second, 

i}) = azimuthal coordinate, positive counterclockwise and 
zero aft, radians. 


n = rotor speed, radians/second. 


Using small angle approximations and the definitions of lift per 
unit length, Z, in terms of bound circulation, r, density, p, and 
tangential velocity, U, 

ji = p u r 

and in terms of the angle of attack, a + V/U, chord, c, lift curve 

slope c and other previously defined quantities, 

Xf cx 

1 - p u2 C c^^ rap + V/U)/2 . 
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By equating these expressions for lift the bound circulation may be 
defined as 


■■ = = “p + ''>/2 


or 


r = c Cjj^^[(ar + sin i|j)ap + V^Ca^ - a^) 


+ w + cos tp + cos sin ^)]/2 


( 7 ) 


where is the sum of collective and cyclic pitch and blade twist, 
radians . 


For a set of rotor blades which are assumed to have stepwise radial 
and azimuthal circulation variations, the above equations may be 
thought of as applying to each radial and azimuthal location. The 
normal component of wake-induced velocity on the blade, w, is made 
up of velocities due to known circulations in the wake and to un- 
known circulations at the blade, and may be v/ritten in the form 

wCr.,4-j^) = Wj,(r. + I I ' '('r) ^ ) (8) 

H 3 

where w (r . ,i|;, ) is the normal component of wake-induced velocity 
^ ^ from known circulations 

r(r.,i|;.) is the blade circulation at (r.,i|/.) 

X, 3 ^ J 

0 .(r.,il^j^) is an influence coefficient which, when multiplied 
^ ^ by the circulation T(r ,i|;.) gives the induced 

K J 

velocity from the (r ,ip.) element of the wake at 

^3 


The summations over the indices £ and j indicate a summation over 
all radial sections of all blades at their respective azimuthal 
positions. Then a set of equations for all r's may be obtained, 
and is of the form 


r(r.,»^) 


c 

2 




£ j 


+ ap(f2r^ + Vj sin + Vf (a^ - a^) + w 


+ ^ cos ip + cos sin ip) 


K 


(9) 
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This equation is solved with a simple iterative procedure (see 
ref. 5 for details) . 


Wake-induced velooity influence coefficients . -Calculation of 
input for the blade loads part of the program is done by adding to- 
gether all a's from the wake program which are multiplied by the 
same bound circulation, r, for computation of the wake-induced 
velocity at a blade load point in the blade loads calculation, and 
by properly subscripting the resulting niimbers. These numbers are 
used in the solution of an equation in the blade loads program 
which is approximately equivalent to equation (9) , except that 
small angles aren't assumed. The evaluation of the input to the 
blade loads program are computed only after a specified number of 
revolutions have been generated. Thus, the input (including both a 
and r type quantities as discussed above) from the wake geometry 
program to the blade loads program is based on approximate (speci- 
fied) blade motions and an approximately repetitive wake. The repet- 
itive characteristic of wake“induced effects is presently determined 
on the basis of previous calculations which depended upon visual 
checking of data and resulting blade loads calculations . 

Large-to-small step size transition . -A considerable savings in 
computer running times and storage was achieved by implementing a 
large-to-small step size option for "final calculations" in the wake 
geometry program. ("Final calculations" are those made during final 
steps, which are taken after the wake is assumed to be repetitive 
and during which wake-induced velocity influence coefficients and 
initial estimates for bound circulations are computed for use in 
the blade loads and response program.) 


Investigation of the forward flight case wake geometries and 
subsequent preliminary calculations indicated that an increase in 
the number of azimuthal steps (or a corresponding decrease in the 
step size) would result in rapidly increasing running times. The 
cost of generating wake program output data for use in maneuver 
loads calculations could be prohibitive for most users. This is 
a consequence of the fact that the number of radial and azimuthal 
blade load points is necessarily much larger for calculations 
which are intended to simulate vortex-blade interactions than 
•those which are not. It was determined during the analysis of the 
forward flight cases that an azimuthal step size on the order of 
10° appeared to be reasonable, and that six or eight radial load 
•points would provide considerably more data than the four used in 
the forward flight case calculations. Also, it was determined 
that blade load response at frequencies through the twelfth har- 
monic of rotor speed be calculated in order to permit comparison 
with measured data, which has response in that frequency range. 

This requirement for calculation of high frequency response neces- 
sitated the use of blade normal modes with correspondingly high 
frequencies. The use of high frequency flapping modes made the 
calculation of air loads at rather closely spaced radial aerodynamic 
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load points imperative in order to compute realistic generalized 
forces. That is, the air load spanwise distribution must have 
variations with the same or smaller spanwise length scale as the 
flapping mode shape variations, at least over the outer portion 
of the blade where high dynamic pressures are expected to occur. 

Preliminary calculations and knowledge of the computer pro- 
gram structure and the types of computer systems which might be 
used to operate these programs led to the following conclusions: 

1. The blade loads program for forward flight could be 
modified for maneuver flight calculations $o as to have up to 8 
radial and 36 azimuthal aerodynamic load points and up to 12 nor- 
mal mode and natural frequency sets without unacceptable increase 
in running time but with a significant increase in the core storage 
required. The increased core storage is a combination of increased 
number of modes and increased number of azimuthal steps per revolu- 
tion. Running time is relatively insensitive to the number of 
modes, and improved programming economics and better modeling of 
control system and lead-lag dampers was expected to improve con- 
vergence and compensate somewhat for the increased number of 
azimuthal steps per revolution. Storage requirements, using over- 
lay methods, was expected to increase from under 70,000 octal to 
approximately 125,000 octal. 

2. The wake geometry program for forward flight could be 
modified for maneuver flight calculations with small increases in 
the amount of core storage required, since only the "full mesh” 
wake had an increased number of elements per blade, and relatively 
few variables ' dimensions were dependent upon the number of azi- 
muthal steps per revolution. However, unless model modifications 
or other changes were made, running times would increase to the 
order of a few hours on the CDC 6600. The running time is a func- 
tion of the number of elements in the wake model, which is approxi- 
mately proportional to the number of trailing vortices (in the 
trailing vortex-only portion of the wake) times the number of 
azimuthal elements in a vortex times the number of blades. The 
running time increases with a power of this number, the power being 
between 3 and 4. The core storage requirements for the blade loads 
and response program was acceptable, but the running time expected ' 
for the wake geometry program was prohibitive for all but research- 
oriented use. As a result of these conclusions, it was determined 
that a change in wake model, or calculation procedure was necessary 

An investigation of wake geometry data for forward flight 
cases indicated that an acceptable approximation to the positions 
of the wake elements of a small step size model (say 36 steps per 
revolution) could be achieved by interpolation between element 
end points of a large step size model (say 12 or 18 steps per 
revolution) . Preliminary calculations confirmed this conclusion . 
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The large-to~small step size option for the wake geometry 
program "final calculations" was implemented in the maneuyer pro- 
gram, according to the following steps, which are indicated in 
figure 4 for clarity. 

1. The program operates as for the forward flight cases 
until (,NxN_"M') steps have been made, and during these calcula- 

tions large step sizes are used, as indicated in figure 4a, Here 

N=l/ii IT, N is the number of azimuthal steps per revolution, and 
s 

N' is defined by 

N' = 1 + of full mesh steps behind each blade-1) 

number of small steps for each large step 

2. The step size is reduced by a factor of 2 or 3, depending 
upon input data, so that by the time the rotor has completed the 

N revolutions the "full mesh" wake immediately behind each blade 
is all small step size, as indicated in figure 4b. All subsequent 
step sizes of the rotor are small. 

3. The "final calculations" are done, with a large portion of 
the trailing-vortex-only wake being large step and a small portion 
immediately adjacent to the full mesh being small, as indicated in 
figure 4c. During these "final calculations" the program inter- 
polates between large step wake element points and effectively 
approximates the large step portion of the trailing-vortex-only 
wake by small step size elements, as indicated by the broken line 
elements in figure 4c. 

All wake program calculations for maneuver cases were done 
with this process. Use of this large-to-small step size option 
is believed to have significantly reduced program running time 
and, in addition, conserved core storage locations. 

Distance-limiting . -The induced velocity from a vortex element 
at a point is approximately proportional to the length of the ele- 
ment divided by the square of the distance between the element and 
the point. Thus the "near neighbors" of a point have the dominant 
effect. Preliminary calculations and investigations of previous 
results indicated that all tip vortex elements further than approxi- 
mately 0.5 to 1.0 rotor radius from a calculation point contribute 
less than 5% of the total wake-induced velocity at a point, yet 
require 50% to 90% of the calculations of this type, for typical 
wakes involved for these rotors and flight conditions. Therefore, 
a distance-limiting option was implemented in the wake geometry 
program. It is felt that this option reduced running time by a 
factor of 2 to 4, and required an insignificant number of additional 
core storage locations. 
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Blade Loads and Response Model and Formulation 

The right-handed coordinate system used in the calculation of 
(periodic) blade circulations and blade response is located such 
that the z-axis is fixed to the shaft, directed upward, the x-axis 
is downstream, blade azimuth angle, ij, is measured with respect to 
the x-axis, y is positive on the advancing side, and the distance 
radially outward from the axis of rotation on a given blade is 
denoted by r , as shown in figure 3 . 

Determination of aerodynamic l oa ding . -The aerodynamic loading 
at a given radial and azimuthal stalilbn is derived from the flow 
experienced by the blade section, as sketched in figure 5* The 
geometric incidence of the section with respect to the rotor plane 
is the sum of the rigid body pitch angle 9^{iJ;) , blade twist r (r) , 

and the torsional deflection <p^(r,jp). The velocity component tan- 
gent to the rotor plane, U, is given by 

U = nr + Vj- cos a sin \jj (10) 

f s 

and that normal to the rotor plane, V, is given by 
V = - V_ sin + h(r,i^) - w ^ w + r(a. cos jp + cos a„ sin ip) 

X S G A K K 

( 11 ) 

where h is the plunging velocity of the section due to the response 
of the blade, and w is the chordwise average wake-induced downwash, 
and is given by 


wtr (,.) = X I o„„(r.,»j) r (r ,♦ ) 

where ^^i ' j ^ wake-induced downwash at (r^,ijjj) due to 

unity bound circulation for blades which were located at 

(r ,4i ), etc. , and w is the Climb rate, 

m n+Ng c 

As is indicated in figure 5, the foil section is replaced by 
a vortex distribution of strength Y(r,<5,ifj) along the chord. This 
distribution is adjusted to make the flow at the section tangent 
to the chord, which relates r to the w and yields the basic rela- 
tionships governing the aerodynamic loading. More details of the 
derivation of r are given in reference 5 . The resulting equation 
for r is 

r(r.,».) = b. (i)^. ,j^2a) 
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where b. = blade semichord, meters, 
= lift coefficient, 
u = (U^ + 


A = 


0^ = 


= 


sum of collective and cyclic pitch angles, radians, 
torsional blade elastic deflection, radians. 


In equation tl2a) the lift coefficient, c„, is a function of angle 

36 

of attack., a, and Mach number, M^, as defined by aerodynamic coeffi- 
cient subroutines. The total angle of attack/ a, is defined by 


cx{r^,il;j) = 6^(ijjj) + (j)(r^,ij;j) + x (rj^) + tan ^ j ) 

(12b) 

where t = built-in blade twist, radians. Equations (12) represent 
a set of nonlinear equations in the strengths', r. The nonlinearity 
is a result of nonlinear aerodynamic coefficient definitions, and 
of the nonlinear dependence of r upon itself (as contained in w and 
its contribution to V, then a, and finally c^) . Solution of equa- 
tions (12) is therefore accomplished in an iterative manner, and is 
discussed in detail in reference 5. 


The lift, drag, and moment per unit span are readily calculated 
once blade circulations have been obtained. Resolving the forces 
into components normal and tangential to the rotor plane, and 

^x' ^respectively , results in the following expressions. 


F = p b • 
ID 


u(U c^ + V c^) 


ID 


Pbi [ 

+ ^ , - r + 2b 


2At j+l D-l 


.(°m ^)j + i " ^°m 


(13) 


P = p b. 
X. . 1 

ID 


u(- V C^ + U c^) 


(14) 


ID 
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where = section drag coefficient, 

c^ = section moment coefficient (about midchord) , 
At = AiJj/n. 

The moment about midchord, Mg, is given by 


Mr 


- pb?|2(<=^ u2) - - r. j). 


1 


3 IT b ? 

1 

8 (At) ^ 


-20 + 0 


’j + i 


9^ 9a_i ^^j + 1 ~ 


(15) 


where 0 

g 


sum of collective and cyclic pitch angles, blade twist, 
and elastic torsional deflection, i.e., the total 
geometric angle of attack, radians. 


5 = local blade spanwise slope, radians. 


The aerodynamic coefficients are computed from empirical approxima- 
tions to experimental data. These empirical approximations are 
given in the Appendix A, for an NACA 0015 airfoil as equations (A-1) , 
in Appendix B for an assumed particular vortex-blade interaction 
model, and in reference 5 for an NACA 0012 airfoil. The aerodynamic 
coefficients are, in general, nonlinear functions of Mach number and 
angle of attack. 


Lumped loads . -Conversion of the aerodynamic loads to a form 
suitable for response calculations is done by lumping the distributed 
loads at the mass points of the lumped parameter blade model. The 
distributed aerodynamic moment is transferred to quarter-chord, then 
the distributed loads are integrated, using a straight-line approxi- 
mation between load points, to obtain lumped loads at the mass points 

The drag force at the i— mass point, for example, is 



^i+l 


F'Cr) dr 

X ^ ^ 


R. 

1 


where R. and R... are midway between masses m. and m. ,, and be- 
X x+ 1 r X— I 

tween m.,, and m. , respectively, and where F^(r) is the straight- 
line approximation to the distributed drag load. Similarly the 
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lumped moment, mg / and lift force f are computed from Mg and F 


z . 

1 


Coordinate transformations are then applied which, result in loads in 
the local blade coordinate system, and are given by 


“ ^z. ^i ■ ^x. ®i 

1. 1. X 


F == f sxn ^ . - f cos 

w. z. X X. X 

XX X 


- irin + z F 
(p , ^ . a . V . 

X X XX 


(16) 


where z = distance of the elastic axis forward of quarter chord, 

^i meters 

$. = angle between chord and rotor plane, positive for airfoil 

^ rotated nose up, radians 

The computation of the total force and moment at each mass, for 
every azimuthal and radial location is done according to the follow- 
ing relationships. 

Q = mQ^ sin $ (h + e cos $) + m9 le + h cos (# - 0g)/cos $] 

V c 

- 2ma’I'e sin $ + + mfg(cos $ cos - sin $ sin sin ip) 

(17 a) 


= mf^^er - 2s50^[Ig sin $ + meh sin 0g/cos $] 
- 2Q4>Iq + 2m0£ (v sin f + w cos $) 


(17b) 


Q = mfi^ cos $ (h + e cos 4>) - mhO sin ($ - 0g)/cos $ 

W G 

- 2maey cos f + F - mf (sin J cos a + cos $ sin a sin t)*) 

w 9^ s s 


(17c) 


Q, = - gj_n $[ig cos 0 + meh] - 0 Ilg + meh cos($ - 0g)/cos $] 

(p c 

+ 2fiij sin f Ig + M - mef (cos $ cos a_ - sin «> sin a sin 


+ f2(a^ cos ij) + sin t|<) Ig 


(17d) 


Qq = n(- sin 4> + cos ij;) (I^ + mr^) 


(17e) 
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where the Q's., i , , v, w and the F ' s and M, are understood to be 

subscripted by i and j ; iti, h, z, r , and I are understood to be 
subscripted by i, and 9 and its time derivatives are understood 

to be subscripted by j. The cyclic pitch, 9^, is given by 

9 = A sin nt + B cos Qt 

c c c 

where = lateral cyclic pitch, radians 

= longitudinal cyclic pitch, radians 

Also, h> is the horizontal separation between the elastic axis and 

the pitch axis at the i— station, and $ . is the total average angle 

between the rotor plane and the chord at the i— radial position, 

th 

and r . is the radius to the i— mass . Reference 5 contains more 
1 

details of equations 17, except that gravity and maneuver load terms 
were not included in the model for the steady level forward flight 
cases. The gravity loading terms are those in equations 17a, c, and 
d which occur as a factor of f^, the gravity loading factor. Maneuver 

rate effects are inclxided in, equations 17d and e, and are those terms 
which are factors of and The maneuver rate effects may be 

thought of as gyroscopic loading terms, with the gyroscopic moment 
being equal to the product of the rotor rotation rate, the appro- 
priate mass moment of inertia, and the precessional (or in this 
instance, the maneuver) angular rate. 

Generalized forces and coordinates . -The generalized forces act- 
ing on each normal mode are computed for each azimuth, according to 




J KJ J 1 


Qv.\. 

X X 


+ Q A 
^W, w. 

X X 


(K) 


'1 11 . 


(18) 


where is from the previous iteration. 


A <>') A 

KJ D • (|) 2 <j) 1 


for J K 


°KK ^ “k 
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d 


K 


is the average aerodynamic damping coefficient for the 
K— mode (read as input) 



is the mode shape quantity for the "q” type of elastic 

deformation, at the i— radial location, for the 
mode 



is the damping associated with motion defined by A 


(K) 


The 0 ^^ terms may be thought of as damping coupling terms. 




The governing equation for the generalized coordinate, 

is given by 




Cr + 


K 


FK(t) 


(19) 


The damping term, is presently defined in terms of the mode 

shape quantities and as 


0 


K 


K 


+ c 



(K) 

4>1 



K 


It should be noted that viscous-type damping of any motion (tor- 
sional, flatwise, or chordwise) may be represented by a similar 
expression. The solution of this equation is obtained in integral 
form. The solution assumes periodicity of both forcing function 
and response. The generalized coordinate, and its time deriva- 

tives are given in Appendix G for an overdamped system, 0 „ > 1, and 

are given in referenGe 5 for an underdamped system, < 1. (It 

should also be noted that the present program is operable for 
0 J, > 1 as well as for 0 ^^ < 1, which was previously not possible.) 

The average aerodynamic damping coefficient term, 29 ( 0 ^, occurs 
on both sides of the governing equation for r (t) . Since an itera- 

XV 

tive solution method is used in computing a compatible set of loads 
and responses, this term effectively cancels at convergence. 


Blade response quantities .- Response variables are computed 

from 

v(t) = I A (t) 

^ V ft. 
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K, 

and similarly for w, W, 41 , and 4 , where 

V = flatwise blade deflection, meters 
w = chordwise blade deflection, meters 

4 > = torsional blade deflection, radians 
= Sv/dr, radians 

V = 3w/3r, radians 

Conversion of the response data to a form suitable for loads 
calculations is then done by computing loads quantities at aero- 
dynamic load positions by linear interpolation of response quanti- 
ties at the mass points. Thus, response quantities, i.e., the plung- 
ing velocity h, the pitch 6, and slope 5 , are given by 

h = Iv cos 4> - w sin $ + z ^ cos 4> + (z - h) (0_ - 05 )] 

3. 3 C 

+ 5V^ cos cos \p , , (20) 

6g = 6o + Oc + <j), and (21) 


5 = - ( 4 cos ■!> + f sin $) . (22) 

These values of h, 6 , and 5 (computed at all aerodynamic radial and 

azimuthal load points) are used in the next iterative calculations of 
blade loads . 

Once convergence has been established between blade loads and 
response, blade shears and moments are computed in terms of the 
normal mode and generalized coordinate quantities according to 

T(t) = I Cj^(t) (23) 

K 

for the torsional moment, and similarly for other moments, shears, 
and motions. 

Determination of initial estimates . -Iterative methods of solu- 
tion are used to solve a nonlinear set of equations, equations (12), 
as well as to determine a compatible set of load and response quan- 
tities. Due to the iterative solution for a compatible set of 
blade loads and blade response, good estimates of load and response 
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quantities are helpful to prevent divergence, and if proper ly chosen 
enhance convergence. A discussion of a particular case Which di- 
verged, and the method used to obtain convergence is discussed in 
the results AND DISCUSSION section. 

Estimates of blade flapping motion for use in blade loads cal- 
culations may be computed by a method equivalent to that used in 
determining blade flapping motions for the wake geometry (see ref. 

5) and wake-induced velocity influence coefficient calculation. 

When the shaft axis is chosen to be perpendicular to the tip path 
plane, blade flapping motions are assumed to be zero. Estimates 
for blade (bound) circulations for the blade loads program are input 
from the wake program, and are those values associated with the last 
steps of the wake program calculations, during which the wake-induced 
velocity influence coeffiGients were computed. 

Blade natural frequency and mode shape model and formulation . - 
The program which calculates blade natural frequencies and mode 
shapes is an independent program which was developed by RASA and 
could be replaced by any approximately equivalent program. The 
necessary output for any such program in order that it might be used 
as input for the blade loads program are the natural frequencies 
together with the corresponding mode shape quantities, i.e., flap- 
ping and edgewise displacements, angles, shear forces, bending 
moments, and torsional deflection and torque. The mode shapes and 
natural frequencies used in the calculations reported here were 
fully coupled edgewise, flatwise, torsion modes. The modes need 
not be fully coupled but should be orthogonal for proper use in 
the blade loads and response program. These mode shape quantities 
are to be defined at the location of the point masses of the lumped 
parameter model, with mode shape magnitudes adjusted to give unity 
generalized mass. The lumped parameter model lengths, masses and 
inertias, mass eccentricities, offsets of elastic axis from pitch 
axis and midchord, twist distribution (as measured at mass points) , 
but not model bending or torsional stiffnesses are used in the blade 
loads-response program. 

Briefly, the model used for the real blade is a lumped parame- 
ter approximation consisting of uniform massless elastic beam sec- 
•tions under tension due to centrifugal loads, with point masses and 
inertias located at the ends of the massless lengths. A modified 
transfer matrix approach is used in determining the natural frequen- 
•cies and mode shapes. (See, for example, ref. 12). The natural 
frequencies and corresponding mode shape quantities are used in the 
calculation of the response of the flexible blades to aerodynamic 
and inertia loads in the blade loads and response part of the pro- 
gram. 
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COMPUTER PROGRAM DESCRIPTION AND USE 


The computer programs used for the implementation of the model 
constitute more than one separate program, and may be classified 
as those used for preliminary ealculations, that used for wake 
geometry and wake influence coefficient calculations, and that used 
for computation of blade loads and blade response. The flow diagram 
shown in figure 1 indicates the general procedure used in the calcu- 
lation of results. The preliminary calculations , some of which were 
done by computer programs, included determination of blade control 
inputs and determination of blade dynamic properties as described 
by a limited n^lmber of fully coupled free vibration mode shapes and 
corresponding natural frequencies. After available flight condi- 
tions, control settings, blade properties, etc. are obtained, this 
data is used to develop input to the wake program and the blade 
loads and response program. Experimentally measured or theoreti- 
cally or empirically calculated data may be used for such data as 
blade cyclic and collective pitch settings, blade flapping motions, 
blade natural frequencies, etc. The blade control settings and 
approximate flapping motions for cases where such data is available 
may be taken from such data. Where experimental data isn't avail- 
able , performance-type programs , charts , or tables may be used to 
provide this data, as appropriate. Blade mode shapes and frequen- 
cies may be computed by any method desired, but must give ortho- 
gonal mode shapes with a generalized mass of unity for use in the 
blade response calculations. The blade vibration characteristics 
used in the calculations discussed herein were obtained with a 
program which computes fully coupled (flapwise, edgewise, torsion) 
free vibration mode shapes and natural frequencies . Airfoil charac- 
teristics used in the calculations were defined in a subroutine by 
algebraic approximations to experimental data. If desired, this 
sxibroutine could be replaced by a similar subroutine or table look- 
up procedure for any airfoil section or for models which include 
rate effects, etc. The calculation or tabulation of airfoil data 
in a form suitable for use in the wake and loads programs may be 
thought of as a preliminary calculation, since at the present time 
it is operational for only NACA 0012 and NACA 0015 airfoil sec- 
tions and one special vortex-blade interaction model , and any 
different models would have to be developed in a similar manner 
before being used. 

After the necessary input quantities have been determined, the 
wake geometry program may be run to obtain wake geometry data and 
wake influence coefficients for use in blade loads. Estimates of 
blade bound circulations are also generated by the wake program for 
use in the blade loads program. 

The blade loads and response program does not require wake 
influence coefficients for operation. Omission results in blade 
loads and response which do not include free wake effects; rather 
an assumed uniform wake-induced inflow is used. It is necessary 
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to input at least one natural frequencY and mode shape from the 
blade free vibration analysis, as well as flight conditions, heli- 
copter blade control settings, and program control variables. The 
blade loads and response program may be conveniently thought of as 
two large sections of which one computes blade loads, and one com- 
putes the response to those loads. Nonlinear effects are included 
in the model, so starting estimates for some quantities are made 
either by input or by the program, and iterations between load and 
response calculations are required in order to compute a compatible 
set of blade loads and response variables. The maximum number of 
iterations allowed is limited by input. Convergence within that 
number of iterations is determined by the relative change in bound 
circulation together with the relative change in motions as mea- 
sured by the generalized coordinates, from one iteration to the 
next, as compared to input numbers. Output of load and response 
data is provided when any of the following situations occurs; 

Cl) the convergence criteria between loads and response is satis- 
fied, (2) divergence of loads from one iteration to the next occurs 
or C3) the maximum number of iterations is reached. This output 
data in the first situation is the desired final output, while in 
the other cases the output serves as a source of information which 
may be used to determine reasons for failure to converge. 


Listings of the program used to compute the wake geometry and 
wake-induced velocity influence coefficients, and the program used 
to compute blade loads and response are given in Volume II, It 
should be noted that the program input and output requires the use 
of English unit quantities, rather than metric unit quantities. 
Metric unit quantities are used in this report. 


RESULTS AND DISCUSSION 


The results of the preliminary calculation are presented and 
discussed in this section, as well as results of the wake program 
and blade loads and response program for maneuver flight conditions 
All calculations of blade loads and response include the effects of 
a free wake. One case includes an empirical definition of lift and 
moment coefficients for a particular vortex-blade interaction model 


Preliminary Ca 1 cu 1 a t i on s 

Preliminary calculations include those necessary to establish 
a numerical model of the helicopter, those necessary to calculate 
blade dynamic properties in terms of normal modes and natural fre- 
quencies, and those necessary to calculate helicopter flight con- 
dition and performance type parameters . 

Helicopter model data . -Information necessary to develop a 
numerical model of the helicopter must be obtained from existing 
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sources. In particular, rotor blade dynamic and aerodynamic proper- 
ties are required. While many rotor blades have the same aerodynamic 
shape, mechanical properties usually differ for each rotor. Suffi- 
cient geometrical and physical properties are needed to permit 
determination of blade dynamic properties , and a more complete dis- 
cussion of normal modes and natural frequencies used in the calcu- 
lations discussed in this report are presented in a subsequent sec- 
tion. When airfoil section data is not in a form suitable for use 
in these computer programs for defining aerodynamic coefficients, 
it must be developed. 

Airfoil section coefficients for the NACA 0015 airfoil were 
developed by using standard curve-fitting techniques to fit low- 
degree polynomiais to experimental data within several sections of 
the angle of attack, range . The source for the experimental data 
was reference 13. Mach number effects were approximated by dividing 
the characteristics as obtained for the low Mach number data of 

reference 13 by the factor ^l-M^j . The resulting equations are 

given in Appendix A. It is suggested that polynomial approximations 
be used as they were here, rather than trigonometric approximations, 
as they are as easily developed, have as great an accuracy, and can 
result in somewhat shorter running times. Rate effects have not 
been modeled for the NACA 0012 or NACA 0015 airfoils, but such 
effects could be included in a similar manner. 


Normal modes and natural frequencies . -It is necessary for oper- 
ation of the blade load$ and response program to have established a 
lumped parameter model of the rotor blades and the corresponding 
fully coupled torsion-f latwise-chordwise , natural frequencies and 
normal modes. (Lumped parameter model details for the H-34 and 
HU-lA helicopters are given in Appendix P.) These calculations 
were performed by an independent program, and resulting normal mode 
and natural frequency data was used as input for the blade loads/ 
response program. Blade natural frequencies and tip mode shape 
quantities and mode identification are given in Tables 1 and 2. 

Mode shapes are given in figures 6, 7, and 8, and include modes with 
frequencies up to approximately 20S^. All of these modes could have 
been used in the blade loads and response calculations. However, 
core storage requirements increase significantly with increased 
number of normal modes. Core storage limitations exist, and present 
a practical restriction on program size. The primary objective of , 
this investigation was the modeling of vortex-blade interactions 
and predicting their effect on aerodynamic loading. Torsional re- 
sponse was expected to be significant for some flight conditions, 
as noted in reference 4. Based upon these considerations, it was 
decided to limit the ntimber of modes. The following criteria were 
used in choosing the modes to be used in the calculations: 

1. torsional response was of considerable interest so "tor- 
sional" modes or modes with significant torsional coupling were 
included. 
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2. "flatwise" modes with frequencies near torsional mode 
frequencies may be coupled aerodynamically and were included for 
this reason, and 

3. modes necessary to deseribe expected steady and low har- 
monic gross blade motions were included, since these were expected 
to have a strong effect on the aerodynamics of the rotor. 

Based on these considerations, modes denoted in Table 1 as F~l, 

F-2, F-3, F-4 , T-1, T-2 , and C-1 were used for H-34 blade repre- 
sentation. Those modes used for representation of the HU-IA teeter- 
ing rotor are those denoted in Table 2 as symmetric modes F-1, F-2, 
F-3 , F-4 , T-1 , and T-2 ; and as antisymmetric modes F-1 , F-2 , F-3 , 
C-1, F-6, and T— 3. The blades operate independently in torsion, for 
the model used in this program, for both rotors. Both teetering and 
fully articulated rotors have pitch bearings and may have torsional 
flexibility because of the pitch control mechanism. In the calcula- 
tions of the H-34 normal modes, the control system stiffness was 
much higher than the blade torsional stiffness, so the torsional 
modes in figure 6 appear to be clamped at the hub. The HU-IA has 
sufficient control flexibility that significant torsional deflec- 
tions occur at the pitch bearing, as may be seen from figures 7 and 
8. Both symmetric and antisymmetric modes are needed to describe 
teetering rotor's flatwise and edgewise response but only pinned- 
free modes are needed for a fully articulated rotor. Thus, more 
modes were required for the HU-IA rotor than for the H-34 rotor. 
Early in the calculations runs were made for the H-34 model with 
modes selected according to the above criteria. Runs were also made 
with more modes in order to evaluate the effect of including addi- 
tional modes. Details of the results of these runs are given in a 
later section. 

Flight conditions and performance parameters .-Parameter values 
such as advance ratio, rotor speed, maneuver rate, air density, 
etc. , were determined from experimental data for those cases where 
experimental data existed. For cases where calculated data was 
used as input, parameter values were chosen which were approximately 
the same as those for similar maneuvers where measured data existed, 
or which seemed to be reasonable values based on nominal helicopter 
maneuver capabilities and flight conditions. In general, the cal- 
culations represent maneuvers which are rather severe, and for some 
of the "calculated" data cases, may exceed helicopter maneuver 
capabilities . 

Blade flap and pitch angles were taken from experimental maneu- 
ver data for cases where available. These values were used in 
calculations for Cases 1 through 6 and 11. It should be noted that 
these measured values did not result in thrust or drag trim. They 
were used because previous calculations indicated that vortex posi- 
tion was a more important factor than vortex strength, and vortex 
positions were expected to be reasonable with use of experimentally 
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measured values of pitch and flap angles. Computed results indicated 
that vortex-blade interactions were predicted, although they were 
not always in agreement with measured data, with respect to phase 
and magnitude. Furthermore, performance programs or charts for 
maneuver flight conditions with nonuniform wake effects were not 
available, and use of the wake and blade loads programs to obtain 
trim conditions would have required significant amounts of computer 
and personnel time with an unknown resulting improvement in pre- 
dicted vortex-blade interaction effects. 

For flight conditions where no experimental data was available, 
approximate values were computed. No performance type computer pro- 
gram or charts were known to be available which readily compute trim 
conditions for a helicopter in steady maneuver flight conditions. 
Computed trim settings were, therefore, determined by using perfor- 
mance charts for forward flight conditions as given in reference 14, 
extrapolating where these charts did not extend into the parameter 
range of a particular maneuver case, and altering these predicted 
forward flight values of pitch settings to approximately account 
for the aerodynamic and gyrosGopic effects of maneuver flight con- 
ditions. Blade flap positions were computed by assuming that the 
shaft was perpendicular to the tip path plane , and that the tip 
path plane was determined by lift and drag trim requirements. 

The calculated and measured rotor parameter values used as in- 
put in the final maneuver case calculations are given in Table 3. 

Of these, the collective angle values for Cases 10 and 12 were re- 
vised after initial blade loads and response calculations to improve 
the predicted thrust values. The blade flapping angles were used in 
the wake geometry program but were not required in the blade loads 
program. All runs were made with the same pitch control settings 
in both the wake geometry program and in the blade loads program. 


Wake Geometry and Induced Velocity Calculations 

The wake geometry numerical computations performed for maneuver 
flight conditions were done by essentially the same procedure as 
those reported in reference 5 for steady state, level, forward flight 
conditions. All cases required two rotor revolutions of wake, and 
running times were approximately 21 minutes CDC 6600 CPU time for 
the four-bladed H-34 and 8 minutes for the two-bladed HU-IA. (Run- 
ning time is strongly dependent on the number of revolutions of 
wake required, as was discussed in more detail in reference 5.) 

All calculations were made with both large-to-small step size and 
distance limiting options, which were necessary to obtain results 
with acceptable run times. 

Wake geometry plots . -Plots of location of the tip vortices of 
the model are useful in visualizing some of the effects of a freely 
distorting wake model. Free wake geometries as predicted by the 
model discussed in this report are significantly different from 
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fixed skewed helical geometries. Resulting wake-induced velocities 
may also differ significantly for a free wake model as opposed to 
a fixed helical wake model, particularly where vortex-blade inter- 
actions occur. A three-dimensional presentation of wake geometries 
would offer the best visualization, but many characteristics may be 
recognized by views such as those presented here . A more detailed 
presentation would include the rear view, and possibly views from 
different angles. It was not the intent of this work to make such 
a presentation. Rather, vortex-blade Interactions are being studied, 
and such effects would not necessarily be more evident in more com- 
plicated or more complete presentations of wake geometry data. It 
is suggested that the most significant visualization for vortex- 
blade interactions would be the presentation of a series of wake 
geometries corresponding to one (or more) rotor revolutions by 
some means such as motion pictures or oscilloscope records. Data 
for such presentations could be most economically and reliably pre- 
pared by’ using computer graphics to display the appropriate wake 
geometry characteristics as determined by the wake geometry program. 
Plots of the wake geometry for an instant in time do not clearly 
display blade— vortex intersections, but they do imply that they 
occur. The use of color could allow more detail or clarity, but isn't 
essential for the purpose of this discussion. These plots do pre- 
sent qualitative information which is useful in the subsequent dis- 
cussion. 

The maneuver wake geometries do not appear to be significantly 
different from level forward flight geometries, especially as seen 
in the top view. Differences in the side and rear views are diffi- 
cult to interpret in terms of vortex-blade interactions . Figure 9 
shows the computed wake geometry for an H-34 coordinated right turn 
at an advance ratio of 0.225, which is referred to as flight 39, 
table 42 in reference 18. Only two tip vortices are shown in the 
figures 9 (a) and (b) and these two plus those portions of the other 
two tip vortices which are near the rotor disk are shown in figure 
9 (c) . The z-scale has been expanded by a factor of five to help 
improve the clarity of the side view plot. It should be noted that 
blade flapping angles are also increased. Inclusion of more vortices 
than these causes unnecessary additional complexity in the plots; 
omission of tip vortices in the rotor disk region would hinder some 
subsequent discussion of this plot. Approximately two rotor revolu- 
tions are presented in figure 9, which corresponds to the time re- 
quired for vortices deposited near \l> = 180° to be transported aft 
of the rotor disk. There is a gross rolling up pattern of the por- 
tions of the tip vortices which were deposited near the sides of the 
rotor disk (see points labelled "a" and "e”) and a humping up of 
those portions of the tip vortices which were deposited near the 
forward portion of the rotor disk (see points labelled "c"). These 
characteristics are very similar to those observed for level forward 
flight cases. 

There are several mechanisms which influence the wake geometry 
pattern, and which combine to give the resulting wake geometry. 
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such as that shown in figure 9. These include the free stream 
velocity, the wake self-induced velocity, induced velocity from 
the bound circulation, and the helicopter maneuver type and rate. 

In figure 9 the direction of the free stream for the rotor is indi- 
cated by the arrow labeled V^. The path of the coordinate system 

origin, which is an indication of the effect of the maneuver, is 
shown by the line labeled "maneuver path". It can be seen that the 
maneuver rate and elapsed time are sufficiently small that the 
maneuver path is approximately a straight line. 

Distortions from skewed helical-type curves in the plan view 
and in the z-direction in the side view are the result of blade- 
wake and wake-wake interactions. The blades appear to come close 
to previously deposited tip vortices on both the retreating and 
advancing sides of the rotor. The tip vortices appear to be inter- 
sected by the rotor blades in both the first and fourth quadrants , 
and are within 0.0 5R over considerable spanwise lengths in the 
fourth quadrant. In the first quadrant the tip vortices appear to 
go from above to below the rotor blades, .{see vortex elements "c-d", 
for example.) In the fourth quadrant the tip vortices appear to 
have one or more local humps and dips. (See vortex elements "a-b-c", 
for example.) As a result, the tip vortex-blade interactions may be 
more severe in the fourth quadrant than in the first quadrant. It 
should be noted that both the first and fourth quadrants have approxi- 
mately parallel orientations of blades and tip vortices, as seen 
from the plan fiew. The primary variables which could result in 
differences of interaction effects are thought to be vortex strength, 
interaction times, separation distance and/or geometrical orienta- 
tion. For normal trim flight conditions the advancing side tip 
vortex strengths will be less than those on the retreating side. 
Interaction times are smaller and frequencies are higher on the 
advancing than on the retreating side due to forward velocity. It 
is felt that these effects, however, are of secondary importance 
with respect to induced velocity magnitude as compared to tip vortex- 
blade separation distance and geometrical orientation. The wake 
geometry shown in figure 10 is for a helicopter in a symmetrical 
pull-up ratio of 0.224. This flight condition was discussed in some 
detail in reference 4, and is referred to herein as Case 6. The 
general characteristics of the two geometries presented in figures 
9 and 10 are similar, i.e., gross roll-up near the edges and local 
humps of vortices which were deposited near = 180® , The free 
stream velocity has an upward component at the rotor disk, which 
is not the normal orientation for steady flight conditions. The 
generally lower position of the downstream vortex elements in figure 
9, as compared to figure 10, is due to the relative change in direc- 
tion of the free stream velocity. This upward component of the free 
stream velocity was present in_ the experimental data for two pull- 
ups and one turn flight condition. Such orientation indicates that 

the helicopter is decelerating in the forward direction, since the 
rotor thrust is perpendicular to the tip path and is acting in the 
drag direction. 
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Wake-induced velocities .-The magnitude and frequency of tip 
vortex-blade interaction effects are more clearly displayed and 
easier to interpret by wake- induced velocity variations tban by 
plots of wake geometry. Wake geometry plots are difficult to inter- 
pret, and more importantly one plot represents only a small fraction 
of th.e time period of one revolution. The use of computer graphies 
could significantly enhance the utility of the wake geometry data 
by displaying the time-varying wake characteristics, as was men- 
tioned previously. (The actual values of wake-induced velocities 
are computed in the blade loads and response program as the pro- 
duct of final bound circulations and wake-induced velocity influ- 
ence coefficients. However, the wake-induced velocity influence 
coefficients used to establish these values are computed by the 
wake program and it is these influence coefficients which determine 
wake— blade interaction. For this reason the calculated wake-induced 
velocity distributions are discussed here rather than in the section 
on blade loads and response results . ) Wake-induced velocities at 
the blade locations are shown for some of the flight conditions in 
figures 11 through 15 for Cases 1, 5, 6, 7, and 10, respectively, 
other flight conditions had wake-induced velocities which were 
generally similar to those in these figures so no others are pre- 
sented. The wake-induced velocities include contributions from all 
wake vortex elements, and are the most distinct representation of 
wake-rotor interactions available. 

Some general characteristics which appear to be consistent for 
all flight conditions are evident from these figures. Maximum wake- 
induced velocities ranged between 0.1 and 0.16 of tip speed, which 
exceeds the inflow component due to V^, for example. The azimuthal 

range for which wake-induced velocities are small and somewhat uni- 
form is generally in the forward portion of the rotor disk and in- 
creases with radius. This is consistent with the fact that with 
forward velocity, the forward portion of the rotor disk has fewer 
vortices near it than the aft portion of the disk, i.e., the blades 
tend to be operating in "undisturbed" air in the forward part of 
the rotor disk. (See figures 9 and 10, for example.) In general, 
the forward third of the rotor disk appears to have very few high 
amplitude interactions in the high dynamic pressure region (out- 
board portion) of the rotor. The first and fourth quadrants show 
significant high frequency, high amplitude variations with some 
large variations occurring in the second and third quadrants near 
90° and 270°, respectively. It should be noted that the peak 
values in these regions may not be the peak values as would be 
predicted by a "continuous" model. Any discrepancies due to the 
"finite time increment" nature of the present model tends to lower 
the peak values, i.e., "peak" velocities may not be computed, but 
may occur between calculation times. It is expected that there 
are not large magnitude variations due to tip vortex-blade inter- 
actions which have a higher frequency than those displayed in 
these plots. Average downstream spacing between tip vortices, 
which is primarily dependent upon forward speed and the rotor 
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speed, are known. These parameters then determine the approximate 
frequency of tip yortex-blade interactions. (See, for example, 
reference 4 , where this frequency is referred to as ) For the 

advance ratios used in these cases the tip yortex-blade interactions 
in the fourth quadrant occur at about the twelfth harmonic of rotor 
speed. In the first quadrant the interaction frequency is higher, 
due to the fact that the vortices and blade are moving in the oppo- 
site directions. Note in figure 9Cc), for example, that there could 
be three interactions in the fourth quadrant and five interactions 
in the first quadrant. (Downstream vortex motion should be taken 
into consideration; vortices move downstream approximately 0.056R 
during one-quarter rotor revolution .) First quadrant interactions 
have wake-induced velocities which have generally lower magnitude 
than fourth quadrant interactions for the same flight condition. 

CThe induced velocity plots shown are machine generated, and are 
straight-line Gonnections of data points which are separated by 
ten-degree azimuth intervals. Thus, frequencies up to the eighteenth 
harmonic of rotor speed are represented.) 

Detailed variations are peculiar to each case, and some of 
these are noteworthy. Cases 5 and 10 (figs. 12 and 15) display 
higher magnitude variations than the others. Case 7 (fig. 14) ap- 
pears to have more interactions in the third and fourth quadrants 
than the other cases. Case 10 (fig. 15) has one strong interaction 
in the fourth quadrant, where a tip vortex was nearly parallel to and 
within approximately 0.03R of the outboard portion (0.5R to l.OR) 
of the blade. This resulted in a wake-induced velocity time history 
which approximates that expected for a parallel interaction with a 
"continuous" time history model of a finite core vortex. This case 
also had the highest computed wake-induced velocity value, 0.1584J2R 
or 34.75 meters/second. 

Effects of wake-induced velocities on total angle of attack . - 
Wake-induced velocities affect aerodynamic loads through their con- 
tribution to the angle of attack. In the calculations, only normal, 
not in-plane, velocity components are computed from wake elements. 
Thus, the total effect of the wake on blade loads is that contained 
in equation (12b), Nonlinear airfoil characteristics, variations of 
the in-plane velocity due to forward velocity, blade cyclic and 
collective pitch, blade flapping, and blade elastic response also 
influence the total angle of attack. Thus, since induced velocities 
are only one of the contributors to angle of attack, some of the 
characteristics of the induced velocity plots are not evident, or 
are suppressed in the resulting plots of angle of attack. This is 
revealed in figures 16 through 19 which show computed angle of 
attack plots for Cases 1, 5, 6, and 10, respectively. Each of these 
figures, when compared with the appropriate one of figures 11, 12, 

13 and 15, illustrate the contribution of wake-induced velocities 
to total blade element angle of attack. (It should be noted that 
positive wake-induced velocities make negative contributions to the 
angle of attack.) Figure 20 shows computed angles of attack for 
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Case 11 (wake-induced velocities were not printed for this case) . 

Wake— induced velocities generally have a reduced effect on 
angle of attack with increased radius, which results from the l/nr 
dependence of arctan (V/U) . High frequency components of angles of 
attack and wake-induced velocities are in phase for all cases. 

Steady and low frequency components of angle of attack are pri- 
marily due to input cyclic and collective. Figures 16, 17, and 
18 are for a four-bladed H-34 helicopter, and display more vortex- 
blade interactions than fig'ures 19 and 20, which are for the two- 
bladed HU-IA helicopter. Figure 19 shows evidence of an interaction 
near ^ = 290 ° with a parallel vortex-blade orientation and small 
distances between the vortex and the blade. No other figures dis- 
play this characteristic . It appears from the plan view of the wake 
geometries that all flight conditions could have near parallel 
vortex-blade orientations, but such interactions are not clearly 
evident in the other figures. The other figures show indications 
of lower magnitude interactions, which resulted from non-parallel 
vortex-blade orientations or larger distances between the vortex 
and the blade than those which were calculated for the flight con- 
dition shown in figure 19. All figures show variations in the 
magnitude or number of interactions with radial location. For 
instance, figure 16 shows more interactions near r - 0.819R than 
at Other radial positions, while figure 17 has one strong inter- 
action at that position and more weaken interactions at other posi- 
tions. The differences for figure 16 are probably due to differing 
distances between the blade and the vortices, with minimum distances 
near r - 0.819R. For figure 17, there appears to have been a vor- 
tex Which went from behind the blade inboard of r = 0 . 8R to ahead 
of the blade outboard of r ~ 0.8R, with a eorresponding phase change 
in the induced Velocity. Figure 17 shows large magnitude inter- 
actions at the inboard positions of the rotor on the retreating 
side and noticeable but smaller and higher frequency variations on 
the advancing side. This figure, when compared with the correspond- 
ing wake-induced velocities of figure 13, clearly displays the in- 
fluence of radius on the contribution of wake-induced velocities to 
angle of attack. Figure 20 shows calculated angles of attack for a 
pull-up flight condition, which appears to have one approximately 
parallel orientation interaction with relatively small magnitude 
effects in the fourth quadrant and several interactions with smaller 
magnitude effects in the first quadrant. 

It should be noted that all calculated angles of attack are 
those determined by an iterative solution of equation (12) , and 
with blade response effects included. Thus, airfoil aerodynamic 
characteristics such as stall have entered into the determination 
of the circulations, and directly affect the wake-induced velocities 
and angles of attack. Different airfoil aerodynamic characteristics 
(for example, effects of angle of attack rates on lift, etc.) would 
result in different bound circulations, different wake-induced 
velocities, and different angles of attack. The angle of attack 
and wake-induced velocity plots indicate that the radial and azi- 
muthal spacing used was sufficient to isolate vortex interactions. 
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to indicate where vortex-induced flow reversed direction, and to 
show expected wake-induced velocity time histories for small separa- 
tion distance near-par allel-type interactions. 


Blade Loads and Response Results 

The blade loads and response calculations performed for the 
maneuver flight conditions were basically similar to those reported 
in reference 5 for steady-state level forward flight conditions, but 
with certain additional considerations. There were gravity loading 
and maneuver gyroscopic loading as included in equations (17) , 
maneuver effects on the velocity normal to the rotor piane as in- 
cluded in equation (11) and, for Case 6, the use of a special tip 
vortex-blade interaction model as previously discussed and as des- 
cribed analytically in Appendix B. Program operation was generally 
the same as for forward flight cases. (A slightly modified calcula- 
tion procedure was used for one case in which a divergence occurred 
in the numerical iteration solution technique related to the solu- 
tion Of equation (12) , which resulted from large wake-induced veloc- 
ity effects. This modified procedure is discussed in the following 
section.) Running times ranged from two to ten minutes (CDC 6600 
CPU time), with most runs taking about five minutes. Run time was 
strongly dependent upon the number of iterations required to obtain 
a converged solution, and was independent of the number of blades. 

The reasons that the run times were independent of the number of 
blades are (1) the program treats only one blade, (2) the same 
number of azimuthal and radial blade aerodynamic load points were 
used for both types of helicopters which were modeled, the H-34 
and the HU-IA, and (3) the number of normal modes used in the re- 
sponse calculation has a negligible effect on the total program 
running time. The data from the wake program which is used as input 
to the blade loads program corresponds to the small-step model dis- 
cussed previously. Radial and azimuthal aerodynamic load point 
spacing was intended to be close enough to display the effects of 
wake-blade interactions, especially as they influence higher har- 
monic loading. Blade load and response quantities which are com- 
puted and printed by this program include aerodynamic loads (forces 
parallel to and normal to the rotor disk, and applied pitching 
moments about the quarterchord) , blade shears, moments, lineal and 
angular response motions , wake-induced velocities (which have been 
discussed previously) and angles of attack. 

Solution for non-convergent case . -The solution of equation (12), 
which is a nonlinear equation for the bound circulations, depends 
upon the accuracy of initial estimates for the bound circulations. 

The normal operating procedure is to use the bound circulation esti- 
mates provided by the wake program as initial estimates for use in 
the blade loads calculation. This procedure was satisfactory for 
the maneuver flight conditions, with the exception of Case 11. This 
case did not converge using this method. The failure to converge 
appeared to be a result of large magnitude wake-induced velocity 
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influence coefficients (which, are indicative of tip yortex-blade 
interactions;) together with a poor Initial estimate for the bound 
circulations. 

The loads program has an option to limit the magnitude of the 
wake-induced velocity effects, so this was done, and a "revised" 
Initial estimate for the bound circulations was obtained by using 
the wake program bound circulations as initial estimates for this 
"revised" set. The limit was then eased somewhat and the process 
repeated, using the "revised" initial estimates to generate a 
"second-level -revised" set of estimates. Although this process was 
continued for several steps , later analysis indicated that a series 
of two revisions of initial estimates (with decreasing limits) would 
have been sufficient to permit convergence. Further discussion of 
Case 11 results are presented in a subsequent section. 

Aerodynamic loads for four— bladed rotor .—Measured section air 
loadings as reported in reference 15 were compared with those com- 
puted. Table 4 gives Case numbers as referred to in this report, and 
flight and table numbers as referred to in reference 15 for the H-34 
helicopter and table and tabulation nxnnbers as referred to in refer- 
ence 16 for the HU-lA helicopter. Both measured and calculated 
blade loads data for some cases for the H-34 helicopter have high 
frequency variations of significant magnitTiide. It is believed that 
these variations may indicate the influence of wake-blade inter- 
actions. Comparing measured and calculated results, it was found 
that the phase and frequency of these variations are generally not 
in good agreement, although the azimuthal region where the inter- 
actions are indicated agree well for some cases. The calculated 
and measured air load variations agree fairly well at r- 0.9 OR for 
the Case 1 flight condition shown in figure 21 , but do not agree as 
well at r-0.75R nor for the flight condition shown in figure 22. 

(In all air load data presentation for the H-34, measured data at 
r=0.90R and r=0. 75R are compared to calculated data at r=0.917R and 
r=0.760R, respectively.) It is clearly evident from the calculated 
angle of attack variations as shown in figures 16 and 17, which 
correspond to these air load variations, that the angle of attack 
variations within approximately ±0.1 radian range are reflected in 
air load variations, while those variations between approximately 
0.1 to 0.3 radians do not appear to affect the air load variations. 
Specific details relating to this effect are discussed later. The 
mean level of the computed section aerodynamic loading is lower than 
the measured loading at all radial positions compared for all cases. 
This overall loading discrepancy is thought to be due to the use of 
measured pitch angles, which resulted in computed thrust values 
which varied between 0.63 and 0.975 of the measured thrust. It 
must be assumed that the helicopter maneuver rates were accurately 
measured. 

The inclusion of Mach number effects on the NACA 0012 airfoil 
characteristics in the static stall angle of attack range greatly 
influenced air load calculations for the H-34 helicopter. 
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Vortex-blade interactions were most common in the first and fourth 
quadrants, as has been discussed previously. The tip Mach numbers 
of the retreating blade were in the range of 0.45 to 0.55, and were 
proportionately smaller with decreasing radius. For this Mach num- 
ber range the lift coefficient is approximately constant for angle 
of attack variations between approximately 0 . 1 and 0 . 35 radians , as 
shown in the sketch below. 



As a result, variations in angle of attack as shown in figures 16 
through 18 have corresponding variations in air load only if they 
range between approximately ±0.1 radian; variations in angle of 
attack between 0.1 to 0 . 35 radians have very small resulting air 
load variations. Retreating blade angle of attack values tend to 
be higher than those for the advancing blade, due to lateral cyclic, 
and were frequently in the 0.1 to 0.3 radian range so that angle of 
attack variations due to wake-induced velocities are generally 
"masked" more on the retreating blade than on the advancing blade. 
While some of the differences between calculated and measured load- 
ings may be due to the low calculated thrust (i.e., incorrect pitch 
angles) , this trend was not observed to be consistent, as seen from 
the comparison of calculated and measured data for other flight con- 
ditions. For example. Case 5, which is a symmetrical pull-up at an 
advance ratio of 0.159 and with a gravity loading factor of 1.6 has 
aerodynamic loads as shown in figure 22. Here, the steady and low 
frequency components appear to be in better agreement than for 
Case 1. It should be noted that Case 5 had a calculated thrust which 
which was 0.975 of that measured- However, at both radial positions 
shown in figure 22 there appear to be more measured vortex-blade 
interactions in the fourth quadrant air loads than were calculated. 
Thus, in terms of the higher frequency aerodynamic loads, it appears 
that Case 1 provides a better representation than Case 5. The 
apparent calculated vortex-blade interactions near t=10°, 110° to 
150°, and 270° to 310° as shown in figure 20 at r=0.75R are not 


35 



indicated in the measured data of reference 15. However, the pres- 
sure time history data at the 9% chordwise station as presented in 
reference 4 for this flight condition indicates that the calculated 
variations at r~0.75H near \})=27Q° did exist. It should be noted 
that the data presented in reference 15 represents a chordwise 
integration of several pressure taps; and that the computed values 
are those determined by ^lifting line theory at the 25% chordwise 
position. The measured data points of reference 15 could lack 
sufficient resolution to show some such high frequency variations 
but it is. felt that other likely reasons for the differences are that 
computed wake positions may be incorrect due to wake program limita- 
tions or pitch and flap angle inaccuracies and that very close 
vortex-blade interactions requires the use of lifting surface theory 
such as that developed (for some special orientations) in reference 
6 . 


In figure 22, there is a general smoothing of calculated air 
loads as compared to measured air loads especially for jjj>210°. The 
angles of attack at these radial positions had significant varia- 
tions in this azimuthal range, as may be seen from figure 17. The 
frequency and phasing of computed angles of attack and measured air 
loads agree much better than do the frequency and phasing of the 
computed and measured air loads , at both radial locations. Angles 
of attack (which have significant contribution froia wake-induced 
velocities) have variations in the fourth quadrant which do not 
influence air loads. The smoothing between calculated angles of 
attack and air loads is primarily due to Mach number static stall 
effects. Similar results were obtained for all other flight con- 
ditions where angle of attack variations occurred within the angle 
of attack range of 0.1 to 0.3 radians and near Mach numbers of 0.5. 
These airfoil characteristic representations effectively "insulate" 
the calculated air loads from the effects of calculated vortex- 
blade interactions. A secondary effect is that the wake-induced 
velocities are determined by using computed bound circulation values, 
and may be somewhat "simoothed" themselves, due to the Mach number 
static stall effects. It should be noted here that no angle of 
attack rate effects were included in the airfoil characteristic 
representations for either the H~34 or the HU-IA helicopter. The 
inclusion of such rate effects could significantly change the 
calculated air loads for cases where angles of attack had relatively 
high magnitudes as well as high frequency variations within the 
static stall regime. NACA 0012 airfoil characteristiG representa- 
tion did not have static stall Mach number effects, as did the H-34 
airfoil characteristics, so there is a one-to-one correspondence be- 
tween angle of attack variation and air load variations for the 
HU-lA helicopter. It appears that such calculated variations may 
generally agree fairly well with measured variations. Therefore, 
it appears that both angle of attack rate effects and Mach number 
effects on stall are important for a realistic representation of an 
airfoil. Furthermore, if rate effects are neglected, omission of 
Mach number effects on static stall may actually improve the repre- 
sentation of airfoil characteristics for vortex-blade interactions. 
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It should be noted that both coordinated turns and syinmetrical pull- 
ups involve a combination of gravity loads greater than one and an 
overall motion which tends to rotate the aft half of the rotor disk 
downward and the forward half upward. 

Measured shaft and flapping angles indicate that both Cases 
1 and 5 are for transient conditions in which the helicopter is 
decelerating in the forward direction. Other calculated and mea- 
sured data were not significantly different than these two cases 
so only a few additional plots of the air loads are presented. 

Calculated and measured air loads for other radial positions 
for Case 1 are shown in figure 23 . There appears to be a computed 
vortex-blade interaction near ^=90° at all radial positions but 
there is decreasing agreement between measured and calculated phase 
and magnitude with decreasing radius. The 9% chordwise pressure 
time history for the radial positions of r=0.55 and outboard all 
have evidence of vortex-blade interactions, and both the phasing 
and the. number of interactions agree fairly well with the calculated 
loadings shown in figure 23- As previously mentioned, the chord- 
wise integration used in reference 15 can result in loss of high 
frequency information which the analysis of reference 4 and the 
calculated data retain. The evaluated high frequency air load 
variations in the fourth quadrant on the outboard portion of the 
rotor are clearly evident in the measured load data of reference 15 
at only the radial positions of r=0.85R and outboard. Again, this 
may be due to vortex position inaccuracies or data integrations as 
done in reference 15, 

Computed and measured air loads for Case 6 are discussed later, 
in the section on use of a particular vortex-blade interaction model, 
and similar results for Cases 10 and 11 are discussed in the section 
on the two-bladed rotor . 


Response of four-bladed rotor . -The response of the rotor as 
measured by the flatwise and chordwise bending moments, the tor- 
sional moment, and the pitch horn load all depend upon applied 
loads (inertial and aerodynamic) and upon blade stiffnesses. Blade 
properties are characterized by the blade normal modes and natural 
frequencies . Where a finite number of normal modes are used to 
represent the rotor blades, blade response is automatically limited 
to that Which can be represented by a linear combination of these 
modes. The time history of the applied loads may or may not be 
properly reflected in the response, depending upon the modes which 
are used to represent the blade. For example, the amplitude of the 
response to a sinusoidal force with frequency, of a particular 

^ — 1 

mode varies with the dynamic amplification factor, [1- , 
where co^, is the K— mode's natural frequency. (The response solu- 
tion described in Appendix C is not derived in terms of forcing 
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function harmonics, but the same comment applies.) Thus, because 
the vortex-blade interaction effects on air loads is usually at 
approxiraatelY the twelfth to fifteenth harmonic of blade passage 
for the H-34, modes with natural frequencies in this frequency 
range could be expected to respond more than others. As was dis- 
cussed previously, torsional response was a primary factor in 
determining the choice of blade modes to be used in the finite set. 
All torsional modes were included, as were modes needed for gross 
blade motions and flatwise modes with significant torsional mass 
coupling or with frequencies close to torsional mode freguencies 
(which could result in significant aerodynamic coupling between 
flatwise and torsional modes) . (See the section on Preliminary 
Calculations for more details.) As a result, flapping modes which 
had frequencies near the twelfth harmonic were not included. This 
resulted in a general suppression of vortex-blade interaction ef- 
fects in the flatwise moments where such variations in air loads 
were computed. 

The flatwise bending moments for Cases 1 and 5 are shown in 
Figure 24, for a radial position at approximately 0.65 of the rotor 
radius. (The measured moment was at r=0.65R and the calculated was 
at r=0,63R for all H-34 flatwise moments.) In general, this is the 
approximate location for the maximum peak-to-peak bending moment for 
an articulated rotor blade. The trend of the calculated flatwise 
moment for Case 1 is in better agreement with the measured data 
than is that for Case 5. The peak values and the phase and fre- 
quency of higher frequency components are not well predicted. The 
suppression of high frequency variations in air loads due to the 
airfoil characteristic representation, as previously discussed, 
effectively eliminated many high frequency air load variations, and 
blade response would not be expected to show the effects of vortex- 
blade interactions in such cases. Case 1, however, had significant 
high frequency variations near ijj=270° on the outboard half of the 
rotor but had no such variations in the computed flatwise bending 
moment, primarily because of the modes used. In general, there were 
significantly more high frequency variations in the measured flat- 
wise moments than in the calculated flatwise moments, The flatwise 
moment values as computed and measured for other cases were gener- 
ally more similar to those presented for Case 1 than for Case 5. 

It should be noted that differences between the physical situation 
and the numerical model could combine to yield Such results. Use 
of a limited set of normal modes and airfoil characteristie repre- 
sentation have been discussed. The relative orientation and loca- 
tion of blades and tip vortices which interact and the timing of 
interactions , as predicted by the wake program is another source of 
differences. Blade air loads and response, which are coupled by 
the nonlinear equations, are also strongly dependent upon pitch 
angles, shaft angle, and other experimentally determined inputs. 

Such input values in the calculations presented here did not yield 
the correct thrusts, and could be expected to also have significant 
effects on other computed values as well. 
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The azimuthal history of computed and measured chordwise bending 
moments at approximately 0.575 of the rotor radius is shown in 
figure 25 for Case 1. (The measured moment was at r==0. 575R, and 
the calculated was at r=0.548R, for the H-34 chordwise moments.) 

The plot shown is typical of measured and calculated chordwise 
moments obtained for all cases for the H-34 rotor, a fully articu- 
lated rotor with a lead-lag damper. The measured data usually 
exhibits a large low frequency (Second or third harmonic of rotor 
speed) oscillation, but the calculated data is nearly uniform and 
positive. Three reasons for this discrepancy are suggested. The 
actual h- 34 lead-lag damper is highly nonlinear. For this damper, 
damping force is essentially proportional to velocity squared up to 
a specified value at which the force is limited by a relief valve. 

The damper constant as contained in reference 15 is considered to 
be too high over much of the damping rate-damping force range. In- 
stead, a damper constant was used which was based on the energy 
dissipated per cycle by the actual nonlinear damper. This is still 
a linear representation, and may be part of the cause of the con- 
sistent discrepancies between measured and calculated values of 
chordwise moment for the H-34 helicopter. The chordwise dynamic 
response of the rotor blades also depends upon the normal modes and 
frequencies used for dynamic representation of the blades. As pre- 
viously discussed, it was desired that high frequency torsional 
response (through the fifteenth harmonic of rotor rotational speed) 
would be permitted. However, core storage requirements precluded 
the use of some of the normal modes computed for the frequency 
range of 0 to 20n. A reasonable criteria was therefore established 
by which some modes were included and other excluded. The same 
criteria were applied to both the H-34 and the HU-lA rotors. One 
of those modes which were excluded was the second chordwise mode, 
denoted C-2, which has a frequency of approximately 3.3n. It now 
appears that this mode had significant response in the measured data. 
While exclusion of this mode may have adversely influenced the chord- 
wise moment calculations , because of the weak coupling involved it 
is not thought to have had significant effect on air load calcula- 
tions or on response displacement-type variables. (Discussion of 
the effects of the use of a larger subset of normal modes is con- 
tinued in the next section.) The other possible reason for the 
large steady chordwise moment is a result of basic limitations of 
the blade frequency data which was used to model the blade flexi- 
bility. The equilibrium position of the blade was assumed to be 
flat and along a radial coordinate in the determination of blade 
normal modes and natural frequencies. The equilibrium position of 
the blade actually involves steady coning and lag. The steady lag 
does not contribute significantly to the chordv/ise moment in the 
real case, but does in the calculations. The assumed equilibrium 
position does not have this large of an effect on the flatwise bend- 
ing moment because the blades are much more flexible in flapping 
than in chordwise bending, and because the coning angle is much 
smaller than the steady lag angle. 
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Torsional moments at mid-span and pitch, horn loads for Cases 1 
and 5 are skown In figure 26- (Both, measured and calculated moments 
are at r=0.5QR, and pitch, horn loads are at the pitch horn for H-34 
data.) The disagreements are significant as considered on the basis 
of general trends, peak values, and high frequency content. The 
peak values as measured In the first and fourth quadrants are under- 
estimated for both the torsional moment and pitch horn load, and 
there Is more high frequency content In the calculated pitch horn 
load within the second quadrant than was measured. These charac- 
teristics were common to all H-34 cases. Both torsional moments and 
pitch horn loads are strongly Influenced by the previously discussed 
airfoil representation Inaccuracies. As may be seen In figure 25, 
there do not appear to be significant frequency .limitations in the 
calculated data. None were expected since the normal modes used 
included all torsional modes. 

Computed and measured torsional moment and pitch horn load data 
are discussed In later sections for the special vortex-blade inter- 
action model and for the HU-IA helicopter. 

Use of additional normal mod es. -More normal modes were calcu- 
lated during preliminary calculatTbns than were used in the blade 
loads and response calculations. The reasons for selection and the 
method of choice of the subset of modes used was discussed in the 
section on preliminary calculations. To investigate the effects of 
including additional modes, one blade loads run was made using Case 6 
wake and loads data, except that additional modes were included. All 
other H-34 simulation runs had been made with seven modes, denoted 
In Table 1 as F-1, F-2 , F-3, F-4, T-1, T-2 , and C-1. The additional 
modes used for the one "Case 6" run only were F-5, F-6, C-2, and C-3. 
The air loads and high frequency response motions appear to have been 
very slightly altered by use of the additional modes. The maximum 
tip peak-to-peak oscillation in the higher modes was approximately 
0.003 meter for flatwise motion and approximately 0.025 meter for 
chordwise motion. The C-2 mode had the largest response of any of 
the four additional modes by an order of magnitude. Because of the 
weak coupling, chordwise motions are the least important with re- 
spect to air loads (with the present model) . The chordwise bending 
moment changed significantly with the inclusion of the additional 
chordwise modes as may be seen in figure 27. It can be seen that 
the frequency and magnitude of the variations of the computed 
moment are in relatively good agreement with the measured varia- 
tions, but that the computed mean value is too high. Except for 
previously discussed items Csuch as possibly incorrect equilibriiam 
position, pitch angles and lead-lag _ damper representation) the 
cause of apparent steady positive bias of the chordwise bending 
moment is not known. No runs with additional modes were made for 
any other cases . 

It is recognized that use of a subset of the normal modes cal- 
culated by a single-blade natural frequency program may adversely 
affect results. Furthermore, the use of normal modes which do not 
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include effects of swashplate flexibility, inter-blade mechanical 
coupling, and other helicopter rotor system characteristics dis- 
cussed in reference 17, for example, may have a significant effect 
on predicted loads and motions. However, within the intent and 
capability of the programs discussed herein, the subsets of modes 
used for the H-34 and the HU-IA appear to allow calculation of the 
dominant air load and response characteristies . (Program modifica- 
tions to include modes such as those determined in reference 17 
would be extensive, since such modes would, in general, be complex, 
and have complex eigenvalues. Extension of the model and program 
described in reference 17 are straightforward, and are being under- 
taken at the time of this writing.) 

Use of special vortex-blade interaction model .- Some of the 
flight conditions for which calculations are reported here were the 
object of an extensive data analysis in reference 4. Based on that 
analysis, empirical relationships for lift and moment coefficients 
were developed for airfoils which undergo stall due to vortex-blade 
interactions. These relationships were implemented as described in 
Appendix B, and were used in computations for one flight condition 
with the blade loads and response program. This case. Case 6 
(referred to as Case 2 in reference 4) , was also computed without 
the special vortex-blade interaction model. Some of the results of 
the calculations are indicated on figures 28 through 30, where three 
sets of data are presented. They are: 

1. measured data as reported in reference 15, 

2. computed data without the special vortex— blade interaction 
effects and with even normal modes used in the blade 
response calculations , and 

3. computed data with the special vortex-blade interaction 
effects and with twelve normal modes used in the blade 
response ealculations . 

The special vortex-blade interactions as modeled used a frequency 
ratio (wj^g/O) of 12, a c^ of 1.0, a Ac^/2 of 0.2, and values of 

c of -0.12. Possible vortex-blade interaction locations 

2D, stalled 

were determined by the magnitude of wake-induced velocity influence 
coefficients, and were restricted to the fourth quadrant. Inter- 
actions were permitted by the VBINT array at the following radial 
and azimuthal locations: for r/R=0.446 at \j)=300° and 350°; for 

r/R=0.597 at i|;=270°, 290°, 300°, 310°, and 350°; for r/R=0.689 at 
i})=270°, 290°, 330°, 340°, and 350°; for r/R=0.760 at t|j= 270°. 290°, 
and 310°; for r/R=0.819 at ij)=270°, 310°, and 320°; for r/R=0.871 at 
t//=280°, 300°, 320°, and 340°; for r/R=0.917 at ^=300°, 310°, and 
350°; and for r/R=0.960 at ip=300°, 310°, 340°, and 350°. The 
interaction criteria, as given in APPENDIX B, were satisfied for 
some of these load positions. The aerodynamic loads as given in 
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figures 28 do show some effects of such interactions. There are 
slightly larger high frequency variations of loading in the first 
quadrant with the special interaction model than without it, but this 
was at least partially due to blade response. (First quadrant air 
loads were not directly influenced by the special model.) There are 
significantly larger high frequency variations in the fourth quadrant, 
which are primarily due to the special vortex interaction. 

Figure 29 shows the computed aerodynamic pitching moments at two 
radial load positions, both with and without the special interaction 
model. The effects of the interaction are clearly evident at 
r/R-0.760, 0.871, and 0.917, but are smaller at 0.960. The data 
presented in reference 4 indicated that strong interactions occurred 
at r/R=0.95. These differences may be due to inaccuracies in the 
wake model calculations or input data, as has been previously dis- 
cussed. Figure 29 does show the frequency and severity of aero- 
dynamic moments which may result from use of the special vortex- 
blade interaction model. Calculations were also made with a value 
of c of -0.24, with larger magnitude but otherwise 

^2D, stalled 

similar air load variations as described above. 

Figure 30 shows the blade response in terms of the blade tor- 
sional moment near r/R=0.50, pitch horn load, and flatwise bending 
moment near r/R=0.55. The vortex-blade interaction effects caused 
small differences between calculated response with and without the 
special vortex-blade interaction model. There are small, high- 
frequency variations in the first and fourth quadrants for the 
special vortex-blade interaction model. The number, location, and 
sign of the variations do not coincide with the measured data as 
given in reference 4. This is due to the locations of possible 
interactions as specified by the VBINT array and other calculation 
model Inaccuracies. The type of response which may be expected from 
the use of the special model is apparent, however. Lack of larger 
magnitude high frequency variations in the torsional moment, as 
compared to the data in reference 4 , may be due to the absence of a 
rotor system torsional natural frequency near the twelfth harmonic 
of rotor speed. Such a frequency could be due to mechanical coupling 
through the swashplate or other inter-blade coupling which was not 
part of the blade model as defined by the normal modes and natural 
frequencies used in this analysis. The flatwise moments shown in 
figure 30 show only low maanitude, high frequency variations with 
the use of the special vortex-blade interaction model. Similar 
characteristics were observed in the data of reference 4. 

The use of the special vortex-blade interaction model appears 
to result in more realistic dynamic stall aerodynamic loads than 
use of static stall characteristics. The use of this or other 
models may be necessary to obtain realistic results for some flight 
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comditions, and requires no significant changes in running time, 
core storage, or input data. Thus, the basic program, and this 
special model, provide a basis for extending the aerodynamic repre- 
sentations to more complex forms than the two-dimensional, static 
stall representation previously used, and indicate the types of 
results which may be expected from such models. 

Aerodynamic loads for two-bladed rotor .- Calculated and measured 
air loads (see ref. 16) and blade response for a symmetrical pull-up 
flight condition are shown in figures 31 through 33. This is the 
only type of maneuver for which comparisons with measured data for 
a two-bladed rotor were made. As with the H-34 rotor cases, the 
rotor control settings for this HU-lA simulation were obtained from 
measured data. The computed thrust was 0.88 of the flight value. 

This thrust value is in better agreement with that measured than 
was obtained for any of the H-34 cases except Case 5. The air load 
levels are' in generally good agreement in the lower harmonics. The 
higher frequency content of the measured data is limited to six 
harmonics , so measured blade-wake interactions are not clearly evi- 
dent, as may be seen from figure 30. It does appear that the calcu- 
lated air loads near 0.75R and 0.90R have tip vortex interactions 
near 90° which the measured data does not contain, and that both 
measured and calculated data indicate an interaction near 300°. 

(The HU-IA calculated air loads near the measured air load positions 
of r=0.75R and r=0. 90R are at r=0.761R and r=0. 919R, respectively.) 

As may be seen in figure 32 (a) , the computed flatwise bending moment 
(measured at r-0. 65R, computed at r-0.63R) has smaller variations 
than those measured, and approximately the same levels; the com- 
puted chordwise bending moment as shown in figure 32 (b) has varia- 
tions with approximately the same magnitude and frequency as those 
measured, but in contrast with the measured data, has essentially 
no steady component. (The chordwise bending moment was measured 
at r=0.60R and calculated at r=0. 59R.) As with the standard H-34 
cases, only one chordwise mode was used in these calculations, that 
denoted the antisymmetric C-1 mode, Which has a frequency of approxi- 
mately 1.5Q. Inclusion of more modes may have significantly changed 
these results. In particular, a symmetrical mode, with frequency 
near zero exists. This would result in an approximately steady bend- 
ing displacement due to steady drag loads, and could account for 
differences in the steady chordwise moment. The torsional moments 
at midspan agree reasonably well in the first and second quadrant 
but the measured results exhibit larger variations in the third and 
fourth quadrant than those calculated, as may be seen in figure 33 (a). 
(Midspan torsional moments for the HU-lA were measured at r=0.50R 
and calculated at r=0,543R.) The pitch loads, as shown in figure 
33 . (b) , have very similar characteristics except that the calculated 
pitch load has a phase lag of approximately 40° compared to that 
measured, and the calculated load generally has a significantly 
higher level than that which was measured. 


43 



Calculated air loads for other HU-IA maneuver conditions were 
obtained, using calculated control angles as previously discussed. 
Both turn conditions had very similar aerodynamic loads and response. 
The air loads, bending moments, and torsional moment and pitch horn 
load are shown in figures 34 through 36 for Case 10. This is a left 
turn at an advance ratio of 0.24 and with a gravity factor of 1.5. 

As may be seen from the air loads near r=0.75R and r=0.90R, the tip 
vortex-blade interaction which was evident in the wake-induced 
velocities (as shown in figure 15) has a more noticeable effect at 
the outboard station. This is consistent with the computed angles 
of attack and airfoil model characteristics. The angle of attack 
has a downward "spike" at 0.9R, but stays below the stall region. 

The angle of attack at r=0.75R goes into the stall region before 
the sudden dip, so air loads show much less effect at that radial 
position than at r“0.9R. 

Elimination of tip vortex .- Rapid dissipation of helicopter 
blade tip vortices has been the object of recent theoretical and 
experimental investigations. (See references 18 and 19.) Elimina- 
tion of the downstream tip vortex in the computer model was done 
for the Case 10 flight condition by defining the tip vortex strength 
in the trailing vortex portion of the model to be zero. Figures 37 
and 38 show computed quantities for this case with and without a tip 
vortex. The severe vortex-blade interaction near ii=290° has been 
eliminated with removal of the tip vortex, with a corresponding 
change in air loads. High frequency air loads in the first quad- 
rant were not eliminated, however. Wake-induced velocities and 
angles of attack in the first quadrant have relatively small varia- 
tions which, together with the rapidly increasing dynamic pressure, 
result in high frequency air load components of significant magni- 
tude. Wake geometry data for the calculation without a tip vortex 
showed that portions of the three inboard trailing vortices which 
were intended to model the vortex sheet have caused the small 
(but effective) wake-induced velocity variations in the first 
quadrant. Blade moments retain some high frequency content, prob- 
ably due to the first quadrant air loads, but are generally 
smoother without than with the tip vortex. 

Vortex blade interaction effect .- The rather strong reductions 
in the vortex-blade interaction effects as observed between the 
computed angles of attack and the computed air loads for the H-34 
cases arises because of the static stall effects as contained in 
the airfoil model. These effects were more severe for some of the 
H-34 cases than for the HU-lA, primarily due to the inclusion of 
more realistic static stall dependence on Mach number for the 
NACA 0012 airfoil than was used for the NACA 0015 airfoil. Dynamic 
effects on the airfoil characteristiGs appear to be necessary in 
order to show the effects of vortex-blade interactions where the 
angle of attack ranges involved extend into the static stall 
region. Significant vortex-blade interaction effects were com- 
puted with the special model used here. 


44 



CONCLUSIONS AND RECOMMENDATIONS 


The wake-induced velocities, angles of attack, aerodynamic 
loads, and blade response have been computed for helicopters in 
steady-state maneuvering flight conditions. Effects of a freely 
distorting wake, flexible blade response, and nonlinear airfoil 
characteristics, were included in the calculations. Computed and 
measured data showed general agreement with respect to the number 
of vortex-blade interactions and their approximate locations, but 
many details did not agree well. The following are specific con- 
clusions obtained from the calculations, from comparison of calcu- 
lated data with measured data, and from the associated analyses 
and interpretations of results. 

1. Tip vortex-blade interactions are predicted in the azi- 
muthal regions, and with frequencies as were expected for heli- 
copters in steady maneuver flight conditions. 


2. The time histories of interactions correspond to that of 
the blade passing near previously deposited tip vortices, which 
occurs primarily in the first and fourth quadrants near the blade 
tips, and over increasing portions of the azimuth with decreasing 
radial position. 

3. Most interactions occur in the first and fourth quadrants. 
First quadrant interactions occur at a higher frequency than fourth 
quadrant interactions, and fourth quadrant interactions commonly 
yield angles of attack which extend into the static stall range. 

4 . Vortex-blade interactions with various orientations and 
distances were computed. Small-distance, near-parallel interactions 
yielded an induced-velocity time history which would be expected 
from a finite core vortex-blade interaction, and occurred primarily 
in the fourth quadrant. Interactions in the first quadrant usually 
had the vortex orientation going from above the rotor disk near 

y=0 to below the rotor with increasing distance starboard.. 

5. Total computed angles of attack clearly displayed wake- 
induced velocity effects, and these effects were approximately 
inversely proportional to blade radial coordinate. Air loads, of 
course, are proportional to velocity squared, so outboard angle of 
attack variations are more important in terms of blade loads than 
are inboard variations. 

6. Calculated aerodynamic loads are dominated by stall 
effects for many conditions where vortex-blade interactions occurred. 
Effects of both angle of attack rate and Mach number must be included 
in aerodynamic representation of the airfoil if high frequency angle 
of attack variations which extend into the stall region are to be 
realistically reflected in computed air loads. 
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7. Due to the fact that reasonably accurate performance 
parameters were not available from performance type programs or 
tables, measured data or data computed by extrapolation and exten- 
sion of level forward flight data was used for the required per- 
formance parameter variables. Measured input data inaccuracies or 
model limitations resulted in low computed thrust values, and yielded 
predicted vortex-blade interactions which seldom agreed in detail 
with measured interactions, as indicated by computed angles of 
attack and measured air loads. 

8. Vortex-blade interaction locations (but not magnitudes) 
show better agreement between computed angles of attack and measured 
pressure tap data than between computed data and measured air loads. 
Lifting surface rather than lifting line representation of the air- 
foil may be necessary for accurate representation of small-distance 
interaction effects. 

It is recommended that further refinements be made and calcu- 
lations be performed which account for the following. 

1. Angle of attack time history and compressibility effects 
should be more accurately represented. Program modifications to 
implement these changes in the blade loads and response would be 
primarily in airfoil aerodynamic representation subroutines. 

2. Lifting surface theory effects should be developed for a 
range of vortex-blade orientations and approximated in so as to be 
able to be efficiently implemented in the existing program. 

3 . Computer graphic representation of wake geometries could 
be used to develop more easily interpreted free wake geometry 
data. Tip vortex-blade interactions could be indicated and time 
histories used to enhance understanding of this complex phenomena. 

4. A performance- type program for calculation of performance 
parameter values for helicopters in steady maneuver flight condi- 
tions should be developed. Such a prograia could include the 
capability to model some of the major effects noted in existing 
studies of helicopter maneuver flight conditions. 
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APPENDIX A 


NACA 0015 AIRFOIL CHARACTERISTICS 


The following equations, referred to as equations (A-1) , are 
based upon NACA 0015 airfoil data reported in reference 3 and 
give aerodynamic coefficients used in the definition of distri- 
buted airloads as a function of angle of attack and Mach number. 

Aerodynamic Coefficients as Polynomial 
Angle of Attack Range Function of Angle of Attack 

0.0 - a - 0.17453 c^^ = 5.4426a 

c^ = 0. 006+1. 2578a2 

c = 1.4467a 
in 

0.17453 < a - 0.59305 c^^ =-5 . 17552006+74 . 39225790a-305 . 08426727a2 

+510 . 05783752a 3-300. 69898135a‘* 

c. = 1.00809827-13. 216159301+60. 24338928a2 

d 

-105. 33915863a 3+66. 65706750a‘^ 

c =-l. 3120201+19. 4541356a-82.038638a2 
m 

+139. 68468 a 3-83. 7789094a'+ 

0.59305 < a - 2.79253 c^ =-l . 1043+5 . 9103a-5 . 2224a2 

+1. 3845a 3-0. 0895a‘^ 

c, =-0. 3690+1. 2002a+1.4200a2 
d 

-1.0989a3+0.1671a‘» 

c =-0 . 1810+1 . 1369a-l . 1730a2 
m 

+0.3019a3-0.0120a'+ 

2.79253 < a - 3.05433 c^^ = 19073 . 5635-26905 . 6199a+14223 . 9143a2 

-3340.04534a3+293.9323a‘» 
c. = 4324. 4195-6164. 2224a+3288.5370a2 

d 

-778. 11679a 3+68. 897280*^ 

c = 20060. 7298-27914. 6308a+14560.552a2 
m 

-3374. 20824a 3+293. 10403a‘^ 

3.05433 < a - 3.141593 c^ =-27 . 057746+8 . 6127481a 

c. = 1. 1263058-0. 110658a2 
d 

c =-13.846155+4. 4073682a 
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For all ranges of a, the lift curve slope, is given by 

dCj/da, and compressibility effects are approximated by dividing 
all quantities by /l-M^ where = Mach number. 
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APPENDIX B 


VORTEX'-BLADE INTERACTION MODEL 


The symmetrical pull-up for the H-34 with advance ratio of 
0.224 was analyzed in reference 4, and empirical relationships 
were developed which approximated the observed lift and moment 
variations. These empirical relationships were implemented accord- 
ing to the following procedure. 


Vortex-blade interactions depend upon tip vortex-blade load 
point proximity and effective angle of attack. Tip vortex-blade 
proximity is reflected in the magnitude of the wake-induced veloc- 
ity influence coefficients , a . It was concluded in reference 4 
that vortex-blade interactions occurred in the fourth quadrant due 
to vortices deposited in the third quadrant. Thus, from these 
experimental observations and by using theoretically calculated o ' s 
an array was established which defined potential vortex-blade inter- 
action locations . If the angle of attack magnitude exceeded a 
specified maximum, an interaction was assumed to have occurred. 

The tests and definitions were done as follows : 


1. Is the current blade load position one with a possible 
vortex-blade interaction (as determined by the input 
array VBINT) ? 


2. Is a > a 7 

— m 

3. If either 1. or 2. are not true, no interaction is indi- 
cated, but if both 1. and 2. are true, an interaction is 
indicated and the lift and moment coefficients are de- 
fined as follows: 


c^(.|.) = + 

u 


AC, 


sxn 


0 ) 


he 


(t - 


(B-l) 
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- 

he 



1 - cos 
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2D, stalled 


(B-2) 


where = angle at which vortex-blade interaction was first 


indicated, i.e. ^ sin ^ 

I 0) 
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u 
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with. c„ lift coefficient calculated from 

normal airfoil cKaracteristics at 

u 


-Tj— = 0.3 Cfrom ref. 4) 


c = two*-dimensional stall limit moment coefficient 

™2d, stalled 



= ratio of vortex interaction frequency to rotor 
speed (referred to as in reference 4) 


c^ =1.0 (for NACA 0012 and NACA 0015 airfoils) 
u 

Once an interaction was indicated, coefficients were defined for 

f 

“he 

one complete cycle by the above method, i.e. , for < 2ir. 
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APPENDIX C 


SOLUTION FOR GENERALIZED COORDINATE OF AN 
OVERDAIIPED SYSTEM 


th 

The governing equation for the K— generalized coordinate, 

“k ~ (G“l) 


is given by 


The solution of this equation is obtained in integral form. The 
solution assumes periodicity of both forcing function and re- 
sponse. The generalized coordinate, and its first time 

derivative, are given below. The second time derivative of 

the generalized coordinate, t"^, is obtained from the governing 

equation itself, after solution for and Subscript K is 

understood on ? , ? , u , ai , 0 , F. 


C(t) = 


e sinh « (t - t) F (x)dT 


Ih 


- 0 “®“^ (Q cosh ut + H sinh cut) 




(C-2) 


k (t) = - 0 UC (t) + 


e x) u)(t - x) F {x)dx 


+ e (G sinh ut + H cosh cut) 


(C-3) 


5 (t) = F Ct) - 2 00 )? (t) - 0)2 ? (t) (C-4) 

where G = (SI + CJ)/(C2 - s^) 

H = (Cl + SJ)/(C2 - s2) 

C = 1 - cosh Z T 

S = e sinh o) T 
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fT _ . > 

I = e cosh &>(T - t) P (T)dT 

- 
Q 

rT 

j = g-aw(T-T) p (x)dT 

0 


(0 


0 ) 



T = 2ir/n 


Computer’ program implementation to evaluate the integrals is done 
by trapezoidal rule. 
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APPENDIX D 


BLADE MODELS DATA 


The following tables give the properties of the blade models 
for the H-34 and HU-IA helicopters which were used to develop 
natural frequencies and normal modes with the blade frequency^ 
program (and which were also used in the blade loads and response 
program calculations) . The tables give the section properties of 
the lumped parameter models, with details of the notation as given 
in references 12 and 13, The lumped parameter blade model begins 
at a radial coordinate at the inboard end of 0.3048 meter (1,0 foot) 
for the H-34 and 0.0 for the HU-lA. The length, bending and tor- 
sional stiffnesses, and mass are those associated with a section 
length- The other quantities are associated with the intersfection, 
or ends, of sections. (The mass is lumped at the ends of sections 
by the computer program, one -ha If of each Section's mass being 
lumped at each end of the section.) 

The bending and torsional stiffnesses are used in the blade 
frequency program but not in the blade loads and response program. 
The Za (or ZA) , which defines the chordwise dimension of the elas- 
tic axis ahead of quarter chord, is used only in the blade loads 
and response program. 


53 



-34 BLADE PROPERTIES 





HU-IA BLADE PROPERTIES 
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APPENDIX E 


SYMBOLS 



B 


b 


C 


C-m 



c 

G 

c 


la 

m 

^ 2 D , stalled 



P 

V 



coordinate transformation matrix (see eq. 2) 
lateral cyclic pitch, radians 

mode shape quantities representing the "q" type 
mode variable (see eq. 18) 

coordinate transformation matrix (see eq. 3) 
longitudinal cyclic pitch, radians 
blade semichord, meters 

generalized coordinate amplitude factor (see 
eq. C-4) 

nomenclature for the m— chordwise normal mode, 
m (an integer) increasing with frequency 

chord, meters 

drag coefficient 

torsional damping coefficient, newton-seconds 
lift coefficient 

lift coefficient limit for special stall model 
(see eq, B-2) 

lift curve slope at 0° angle of attack 
moment coefficient 

two-dimensional stall limit moment coefficient 

coordinate transformation matrix (see eq. 2) 

generalized force for the K— mode, newton- 
meters 

lumped aerodynamic forces acting at blade mass 
points, newtons 
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distributed aerodynamic forces per unit span, 
newtons/meter 

+* 

nomenclature for the m— flatwise normal mode, 
m (an integer) increasing with frequency 

gravity load factor, thrust/weight 

generalized coordinate amplitude factor (see 
eq. C-4) 

helicopter acceleration, meters/second^ 

generalized coordinate amplitude factor (see 
eq. C-4) 

chordwise dimension of the elastic axis ahead of 
the pitch axis , meters 

blade response plunging velocity, meters/second 

generalized coordinate amplitude factor (see 
eq. C-4) 

blade element torsional mass moment of inertia 
about the elastic axis and c.g., respectively, 
kilogram meters^ 

blade element chordwise mass moment of inertia 
about the c.g. (I hi), kilogram meters ^ 

generalized coordinate amplitude factor (see 
eq. C-4) 

lift per unit span, newtons/meter 
Mach number 

distributed aerodynamic moment per unit span, 
newtons 

lumped aerodynamic moment acting at blade mass 
point, newton/meters 

blade element mass, kilograms 

number of revolutions, N - l/yir 

number of azimuthal steps per revolution 


coordinate transformation matrix (see eq. 2) 



q 

R . 


Rc-Rp 


r 

S 

T 


T-m 


t 

U 


u 

V 



w 


x,y,z 




m 




total "q-type" load at blade element mass point 

coordinate transformation matrix (see eq. 1) , or 
rotor tip radius, meters 

maneuver radii for coordinated turns and sym- 
metrical pull-ups, respectively, meters 

radial coordinate, meters (see fig. 4) 

coordinate transformation matrix (see eq. 1) or 
generalized coordinate amplitude factor (see eq.c-4) 

torque, newton/meters (see eq. 22) ; or period, 
seconds (see eq. C-4) 

th 

nomenclature for the m— torsional normal mode , 
m (an integer) increasing with frequency 

time, seconds 

velocity tangential to rotor plane, meters/second 
(see eq. 6) 

(U^+V^)^'^^, meters/second 

velocity normal to rotor plane, meters/second 
(see eq. 5) 

helicopter velocity, meters/second (see fig. 4) 
climb rate, meters/second 

wake-induced velocity normal to rotor plane, 
meters/second 

rectangular coordinates, meters (see fig. 4); 

X is positive aft, 

y is positive on the advancing side, and 
z is parallel to the shaft and positive up 

total blade section angle of attack, radians 

vortex interaction angle of attack limit, radians 

pitch, angle increment, radians Csee eq. 4) 

stall angle of attack (see APPENDIX B) , radians 

sum of pitch and blade twist angles, radians 

angle between the roll axis and the x-axis, 
radians 
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Bo 

Bi 


32 

r 

Y 

6 

* 

6 

e 


e 

g 


9o 


0 


T 


41 

A 

A 


y 


P 

a 


shaft tilt angle, positive aft, radians 

longitudinal first harmonic blade flapping angle, 
positive for flap up at \jj=0 , radians 

steady blade coning angle , radians 

longitudinal first harmonic flapping angle, 
radians 

lateral first harmonic flapping angle, radians 

circulation, meters /second 

circulation per unit chord, meters/second 
(see fig. 5) 

chordwise coordinate, meters (see fig. 5) 
turn rate, radians/second 

distance of mass c.g. ahead of elastic axis, 
meters 

th. 

generalized coordinate of the mode 
cyclic pitch, radians 

total geometric angle of attack, 0^ + x + (|)^, 
radians 

collective pitch, radians 
So + 0c 

3v/8r, radians 

incremental turn angle, 6 At, radians (see eg. 2) 
bound circulation term (see eg. 12) 

advance ratio, V^/nR 

coordinate transformation matrix (see eg. 3) 

blade spanwise slope, radians/meter 

air mass density, kilograms/meter ^ 

wake-induced velocity influence coefficient, 
1/meter (see eg. 8) or damping coefficient 
(see eg. 18) 
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a 

T 

$ 


'f 

'P 


he 


Subscripts : 


A 

i 

J 

j 

K 

k 

£ 

n 

0 

q 

R 

V 


average aerodynamic damping coefficient 
blade twist, radians 

angle between rotor plane and chord, positive 
for airfoil nose-up, radians 

coning angle for turns, radians (see eq. 2) 

incremental roll angle, radians (see eq. 3) 

torisonal elastic deflection, radians 

chordwise elastic bending slope, radians/meter 

azimuthal angle, radians (see fig. 4) 

rotor speed, radians/second 

blade natural frequeny 

vortex interaction frequency 

w/l- 0 ^ , or co/o^-l for underdamped or over- 
damped systems, respectively, radians/second 
(see eq. C-4) 

denotes pitch maneuver quantities 

denotes radial location 

normal mode number 

denotes azimuthal location 

normal mode number 

denotes azimuthal location 

denotes radial location 

denotes radial location 

denotes a steady condition or parameter 

a dummy subscript variable (see eq. 18) 

denotes roll maneuver quantities 

denotes flatwise blade load and response quan- 
tities 
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denotes chordwise blade load and response quanti- 
ties 

denotes quantities parallel to the x-axis or 
normal to the rotor shaft 

denotes quantities parallel to the z-axis 

denotes spanwise slope type blade load and 
response quantities 

denotes torsional blade load and response quan- 
tities 

denotes chordwise blade load and response quan- 
tities 


Mathematical Notation 

denotes an increment, i.e.. At is an increment 
in t 

indicates time derivative, i.e., v=dv/dt 
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TABLE 1 


H-34 FREQUENCIES AND MODE SHAPES 
n=23.0 rad/sec, 00=13.6° 
FULLY ARTICULATED 


0) 

Tip Def; 

Lection 

Mode 

Coupling, 
if important 

rad/ sec 

X 1/fi 

<j>/v 

-w/v 

Type 

5.6804 

0.2470 

-0.0110 

-9.60000 

C-1 

F-1 

23.6361 

1.0277 

-0.0082 

0.09230 

F-1 

0 

61.2762 

2,6642 

-0.0010 

-0.11300 

F-2 

0 

75.6725 

3,2901 

-0.0163 

5.39000 

C-2 

F-3 

111.3324 

4.8405 

0,0114 

-0.08870 

F-3 

0 

163.2693 

7,0987 

780.8 

0.63900 

T-1 

0 

171.3437 

7.4497 

0.1345 

0. 12400 : 

F-4 

C-3 

192.8873 

8.3864 

0.1175 

-8,05000 

C-3 

F-3 

250.5688 

10,8943 

-0.0301 

0.00925 

F-5 

0 

347.3836 

15.1036 

0.0224 

0.05620 

F-6 

0 

357.5299 

15.5448 

0.0129 

-15,9300 

C-4 

0 

444.9662 

19.3464 

-0.0928 

-0,01510 

F-7 

0 

484.3651 

21,0594 

57.74 

-0.01160 

T-2 

0 
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TABLE 2 


HU-IA FREQUENCIES AND MODE SHAPES 
fl=32.882 rad/sec, 0o=lO° 


0) 

Tip Deflection 

Mode 

Coupling , 

rad/sec 

X 1/n 

<p/v 

-w/v 

Type 

if important 

symn 

etric mode 

s (clamped 

, flap, pin 

ned choi 

rdwise) 

37.1260 

1.1291 

-0.0994 

-.06200 

F-1 

0 

104.9312 

3,1910 

-0.8530 

-.06330 

F-2 

T-1 

126,5499 

3.8486 

15. 0600 

-.06080 

T-1 

F-2 

181.6353 

5.5239 

-0.8464 

2.89000 

C-2 

F-3, T-1 

200.2473 

6.0899 

-0.9613 

-.34400 

F-3 

T-1, C-2 

310.8473 

9.4534 

1.8440 

+.01520 

F-4 

T-2 

338.2845 

10.2878 

15.6800 

-.01300 

T-2 

F-4 

468.2712 

14.2410 

0.6758 

-.05780 

F-5 

T-3 

546.1555 

16.6096 

-0.1732 

11.27000 

C-3 

0 

621.0826 

18.8882 

-13.470 

-.01220 

T-3 

F-6 

636.4494 

19.3556 

-1.1700 

-.06120 

F-6 

T-3 

antisyrametric inodes (pinned flap, clamped chordwise) 

32.0184 

0.9737 

-0.0725 

-.34400 

F-1 

T-1 

4 9.7393 

1.5127 

-0.1883 

2.80000 

C-1 

F-1 

90.3910 

2.7490 

-0.6168 

-.19000 

F-2 

T-1, C-1 

126.4704 

3.8462 

21.6900 

-,07820 

T-~l 

0 

161.8563 

4.9223 

-0.9383 

0.01600 

F-3 

T-1 

271.1910 

8.2474 

0.0892 

0.19700 

F-4 

c-2 

288.7094 

8.7802 

0.2464 

-4.39000 

C-2 

F-4 

339.2609 

10.3175 

19.7800 

-.05150 

T-2 

F-4 

414.4594 

12.6044 

-0.3169 

+.03120 

F-5 

T-2 

601.8897 

18.3045 

4.081 

-.01107 

F-6 

T-3 

609.0193 

18.5214 

6.378 

-.01620 

T-3 

F-6 
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ROTOR PARAMETER INPUT VALUES 
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Cases 1 through 8 are for an H-34 heliGopter, Cases 9 through 12 are for an HU-lA. 
q is a for pull-ups and turns, for rolls. 




CASE NUMBERS AND DATA SOURCES 
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Required thrust is that measured, where reported; and that which would result 
from the specified load factor for "computed" cases. 





Figure 1. Program usage flow diagram 




Figure 2. Wake model with combination of "full mesh" 
wake of shed and trailing vortices and 
"modified" wake of trailing vortices only 
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Figure 3, Rotor system coordinates 
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(a) "large" steps only 



(b) wake repetitive/ ready 
for "final calculations, 
with full mesh wake made 
up of "small" steps 


ations" 


Large-to-small step size procedure, illustrated here 
for eight large steps and sixteen small steps per 
revolution. All calculations were performed with 
eighteen large steps and thirty ^six small steps per 
revolution. 




Flatwise Bending Modes Chordwise Bending Modes 

Type Type 

m/n flatwise deflection m/n chordwise deflection 


F-1 

1.03 


F-2 

2.66 


F-3 

4.84 


F-4 

7.45 


F-5 

10.9 


F-6 

15.1 


F-7 

19.3 



Figure 6. H-34 coupled mode shapes and frequencies, 

n * 23 rad/sec (Very weak coupling exists 
between mode types, so only major mode shape 
is shown and all mode shapes are for a tip 
flatwise deflection of unity.) 
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Figure 11. 


Computed wake-induced velocities at 
aerodyncunic load points for Case 1, 
right turn, p=0.185, f =1.5 
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aerodynamic load points for Case 10, 
HU-IA, left turn/ y=0.24, f^=1.5 
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computed angles of attack, radians 






degrees 


Figure 17. Computed angles of attack at aerodynamic load 
points for Case 5, H-34, pull-up, u==0.159, f = 
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Figure 20. Computed angles of attack at aerodynamic load 
points for Case 11, HU-IA, pull-up, y-0.21, 
f =1.15 
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Figure 22. Section aerodynamic loading for Case 5 
pull-up, u=0.159, f =1.6 
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chordwise bendxnq moment 



pitch horn load torsional moment pitch horn load torsional 
newtons newton meters newtons newton : 


O measured, ref 
— computed 


^ , degrees 

(a) Case 1, torsional moment 


w, degrees 

(b) Case 1/ pitch horn load 


(c) Case 5, torsional moment 




ip , degrees 


(d) Case 5, pitch horn load 


Fxgure 26. Torsional moments at 0.5R and pitch horn 
loads for Cases 1 and 5 
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(b) loading for r=0.90R 

Figure 34. Section aerodynamic loading for Case 10 
HU-IA, left turn, y=0.24, f =1.5, 
computed only 
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Figure 37. Case 10 convputed wake-induced velocities and 

angles of attack with and without tip vortices 
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